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Electron  Distribution  in  Molecules.  IV.  Phosphorus  Magnetic  Resonance  Shifts* 

H.  S.  Gutowsky  and  D.  W.  McCall! 

Noyes  Chemical  Laboratory,  University  of  Illinois,  Urbana,  Illinois 
(Received  August  28,  1953) 

Chemical  shifts  in  the  phosphorus  nuclear  magnetic  resonance  absorption  have  been  measured  in  twenty 
compounds.  The  shifts  arise  from  differences  in  the  magnetic  shielding  of  the  phosphorus  nuclei  by  elec¬ 
trons,  and  the  values  observed  are  interpreted  qualitatively  in  terms  of  the  molecular  electronic  structures. 

The  range  in  shifts  is  7X10-4  which  is  somewhat  larger  than  observed  for  fluorine  on  which  more  extensive 
data  are  available.  The  phosphorus  resonance  shifts  fall  into  two  distinct  groups  one  of  which  contains 
compounds  with  trivalent  phosphorus  and  the  other,  pentavalent.  The  shifts  in  pentavalent  compounds  fall 
in  a  relatively  narrow  range,  generally  with  greater  nuclear  magnetic  shielding  than  in  the  trivalent.  The 
resonance  shifts  between  similar  or  homologous  compounds  can  be  accounted  for  by  considering  changes  in 
the  importance  of  ionic  and  double-bonded  electronic  structures. 


INTRODUCTION 

VARIATIONS  in  the  location  of  a  nuclear  magnetic 
resonance  with  chemical  constitution  were  first 
reported  by  Knight1  for  P31.  Although  no  details  were 
given,  his  estimated  magnitude  of  0.01  percent  for  the 
effects  is  consistent  with  the  measurements  reported 
herein.  Ramsey’s2  theoretical  analysis  ascribes  the  varia¬ 
tions  to  magnetic  shielding  of  nuclei  by  the  molecular 
electrons,  whose  motions  interact  with  the  applied 
magnetic  field.  A  change  in  electron  distribution  with 
chemical  constitution  changes  the  magnetic  shielding 
of  a  nucleus  and  chemically  shifts  its  resonance  position. 
Thus,  in  principle,  information  about  the  electronic 
structures  of  molecules  can  be  obtained  from  chemical 
shift  data.  Although  detailed  calculations  have  been 
made2  only  for  H2,  data  and  semi-empirical  analyses 
have  been  given  for  a  large  number  of  proton  and 
fluorine  chemical  shifts  in  various  types  of  compounds.3-6 

If  the  chemical  shifts  are  to  be  measured  accurately, 
the  resonance  lines  should  be  narrower  than  the  shifts, 
and  this  usually  requires  that  the  sample  be  liquid  to 
minimize  internuclear  magnetic  dipole  broadening. 
Also,  as  the  resonances  can  be  broadened  by  electric 
quadrupole  interactions,  the  best  nuclei  for  this  type 
of  experiment  are  those  with  spins  of  §,  in  which  case 
the  electric  quadrupole  moment  is  zero.  Furthermore, 
the  nuclei  must  exist  in  sufficient  isotopic  abundance 
and  with  a  large  enough  magnetic  moment  to  give 
measurable  absorption.  Fluorine  and  hydrogen  nuclei 
are  well  suited  for  this  type  of  experiment  as  evidenced 
by  the  results  already  obtained.  Phosphorus  has  one 
stable  isotope  PM,  and  it  meets  the  requirements  for 

*  Supported  by  the  U.  S.  Office  of  Naval  Research, 
t  National  Science  Foundation  Predoctoral  Fellow. 
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shift  measurements  rather  well  as  it  has  a  spin  of  and 
a  sizeable  magnetic  moment. 

The  main  objective  of  this  note  is  to  summarize  the 
phosphorus  resonance  shifts  observed  thus  far  and  to 
discuss  their  utility  in  determining  molecular  differ¬ 
ences  in  bond  type.  Some  of  the  phosphorus  resonance 
shifts  we  have  observed  were  given  in  connection  with 
the  multiplet  resonance  lines  arising  from  the  electron 
coupling  of  nuclear  spins.7  Also,  Dickinson8  reported 
earlier  the  resonance  shifts  in  PCh,  PBr3,  and  POCl3 
with  respect  to  H3PO4  (aq.) ;  our  results  agree  with  his 
within  experimental  error.  In  considering  the  phos¬ 
phorus  resonance  shifts  we  encounter  two  factors  not 
present  in  case  of  fluorine  and  hydrogen.  First,  phos¬ 
phorus  is  multivalent  and  second,  it  has  two  stable 
valence  states.  Although  these  factors  complicate  any 
interpretation,  one  may  be  led  to  a  better  understanding 
of  multivalent  bond  structures  in  the  attempt. 

EXPERIMENTAL  METHOD 

The  experimental  details  were  similar  to  these  de¬ 
scribed  earlier.3’7  The  resonances  were  observed  at 
about  6365  gauss,  using  an  rf  bridge,  30  cps  modulation, 
and  oscilloscope  display.  The  resonance  shifts  were 
measured  at  a  fixed  radio-frequency  as  the  difference  in 
applied  field  to  center  the  resonance  lines  of  the  com¬ 
pound  and  the  reference  on  the  oscilloscope.  The  multi¬ 
plet  structures  observed  in  many  of  the  resonances  were 
symmetrical  and  the  center  was  used  in  determining 
the  resonance  shift,  as  is  justified  theoretically.7  Usually 
a  set  of  ten  measurements  was  made  by  placing  suc¬ 
cessively  the  sample  and  reference  compound  in  the 
same  coil.  The  order  was  reversed  for  half  the  measure¬ 
ments  so  any  constant  drift  in  frequency  or  field  would 
cancel  out.  Statistical  errors  were  the  order  of  ±0.006 
gauss,  and  systematic  errors  up  to  one  percent  may  be 
present  in  addition.3 

Compounds  do  not  need  to  be  extremely  pure 
for  chemical  shift  measurements,  and  commercial 

7  Gutowskv,  McCall,  and  Slichter,  J.  Chem.  Phys,  21,  279 
(1953). 

S\V.  C.  Dickinson,  Phys.  Rev.  81,  717  (1951). 
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samples  were  used  without  further  purification.  The 
C1(CH,);iPC12  and  C2H5OPOCI2  were  kindly  provided 
by  Dr.  D.  R.  Martin;  the  sources  of  the  other  samples 
have  been  given.7  The  shifts  in  PH3  and  PF3  were 
measured  in  the  liquids  at  about  —  90°C.  Yellow  phos¬ 
phorus  (P<)  and  P4S3  were  measured  at  room  tempera¬ 
ture  in  the  solid  state,  the  latter  with  a  narrow  band 
amplifier.  The  other  compounds  were  liquid  at  room 
temperature  and  were  measured  as  such. 

In  this  paper  we  define  the  chemical  shifts  of  P31  as 
105X  (Hc—Hr)/Hr  where  Hc  and  II r  are  the  magnetic 
fields  required  for  resonance  in  the  sample  and  in  a 
reference,  respectively,  at  a  fixed  radio-frequency.  Aque¬ 
ous  H3PO4  was  chosen  as  an  arbitrary  reference,  as  it 
gives  a  strong,  narrow  resonance  and  is  readily  available. 
With  this  definition,  more  positive  shifts  correspond 
to  greater  magnetic  shielding  of  the  nucleus.  The 
3  values  used  in  discussions  of  fluorine4'6  and  proton5,6 
shifts  are  defined  similarly,  but  the  sign  is  reversed  so 
that  changes  in  S  correspond  to  changes  in  the  charge  on 
the  atom. 


RESULTS  AND  DISCUSSION 

The  phosphorus  resonance  shifts  observed  are  given 
in  Table  I,  with  trivalent  and  pentavalent  compounds 
listed  separately.  Before  discussing  individually  several 
of  the  various  types  of  compounds,  some  general  ob¬ 
servations  will  be  given. 

The  pentavalent  compounds,  excluding  HPF6,  have 
phosphorus  shifts  clustered  in  the  narrow  range  from 
2  to  —3.  In  HPF«,  which  has  a  higher  shielding  of 
11.8,  the  structure  is  octahedral,  spW,  while  the  other 
pentavalent  compounds  use  sp*d  orbitals.  In  Ramsey’s 
theory2  the  shielding  increases  with  the  electrostatic 
potential  at  the  nucleus  and  is  reduced  by  a  second- 
order  paramagnetism  which  arises  in  part  from  the  lack 
of  spherical  symmetry  introduced  by  the  bonding  elec¬ 
trons.  The  higher  symmetry  of  the  PF6~  compared  to 
the  other  compounds,  as  well  as  the  contribution 
of  the  additional  two  electrons  may  account  for  the 
greater  shielding  of  the  phosphorus.  Also,  as  we  shall 
consider  below,  replacement  of  other  atoms  with  fluorine 
increases  the  phosphorus  shielding,  and  PF#-  is  unusual 
in  the  number  of  fluorines. 

With  the  exception  of  P4,  PH3,  and  P4S3,  the  trivalent 
phosphorus  compounds  have  shielding  values  ranging 
from  —9.7  to  —22.2,  considerably  below  the  —3  lower 
limit  of  the  pentavalent  compounds.  The  bonding  in  the 
trivalent  compounds  is  a  combination  of  p 3  and  sp~, 
and  the  lower  shielding  most  likely  arises  because  there 
are  fewer  valence  electrons  and  a  less  symmetrical 
electronic  distribution  than  in  the  pentavalent  com¬ 
pounds.  In  this  connection  it  is  of  interest  to  note  that 
the  chlorine  nuclear  quadrupole  coupling  constants2'"' 


Tablk  I.  Chemical  shifts8  in  the  P31  magnetic  resonance. 


Trivalent  compounds  IVntavalent  compounds 


Formula 

Chemical 

shift 

Probable 

error 

Formula 

Chemical 

shift 

Probable 

error 

P,  (in  CSol 

48.8 

0.02 

HPFo 

11.8"'' 

P, 

45.0 

0.03 

F2PO(OH) 

2.01" 

0.01 

PHs 

24.  1" 

0.03 

POCIFa 

1.48" 

0.02 

PlSs 

11.4° 

POCljF 

0.00" 

0.02 

PFs 

—  9.7" 

0.05 

H3PO4 

0.00 

reference 

CHjOPF. 

-11.1" 

0.02 

HPO(OH). 

-0.45" 

0.04 

PaL  (in  CS2) 

-17.0' 

POCls 

-0.54 

0.04 

Pis 

-17.8 

0.04 

C2H50POC12 

-0.64 

0.06 

C1(CH,)3PC1: 

>  -18.2 

0.07 

HjPO(OH) 

-1.38" 

0.06 

PC13 

-21.5 

0.09 

PSCh 

-3.08 

0.05 

PBr3 

-22.2 

0.11 

n  Referred  to  aqueous  HjPCh;  defined  as  10 ’•X(He-'Hr)/Hr  where  lie 
and  Hr  are  the  resonance  magnetic  fields  for  the  sample  compound  and  refer¬ 
ence,  respectively. 

b  The  resonance  exhibited  multiple!  structure;  see  reference  7  for  details. 
0  The  resonance  was  weak  and  the  shift  value  is  good  only  to  ±1.0. 


demonstrate  that  the  P—  Cl  bonds  are  considerably 
more  covalent  in  PCI5  than  in  PC13. 

(a)  PF3,  PC13,  PBr3,  PI3 :  In  this  homologous  series 
the  phosphorus  shielding  decreases  in  the  order  PF3, 
PI3,  PC13,  and  PBr3,  with  the  iodide  out  of  sequence 
and  the  fluoride  shielded  considerably  more  than  the 
other  halides.  Similar  trends  have  been  noted  in  the 
proton  shifts6  in  the  methylene  halides  (CH2Y2)  and 
haloforms  (CHA'3),  and  in  the  substituent  effects4  in 
benzene.  The  observed  order  reflects  the  presence  of 
two  opposing  effects,  ionic  and  double  bond  character. 
The  P— X  electronegativity  difference  increases  on 
passing  from  PI3  to  PF3,  and  the  ionic  bond  character 
should  increase  accordingly.  Also,  fluorine  and  to  a 
lesser  extent  chlorine  and  the  larger  halogens  can  form 
double  bonds.  So  the  observed  shifts  are  compatible 
with  the  phosphorus  shielding  increasing  as  the  double 
bond  character  is  increased  and  decreasing  as  the  ionic 
bond  character  is  increased.  Thus,  for  ionic  character 
effects  alone,  the  phosphorus  shielding  would  decrease 
in  the  sequence  I,  Br,  Cl,  and  F ;  the  opposing  double 
bond  effect  is  greatest  for  F  and  the  next  for  Cl,  giv¬ 
ing  the  observed  sequence  F,  I,  Cl,  and  Br.  More¬ 
over,  the  double  bond  character  effect  in  PF3  is  very 
pronounced  compared  with  the  other  halogens;  this 
is  consistent  with  the  well-known  chemical  concept 
that  first  row  atoms  double  bond  most  effectively." 

(b)  P0C13,  POCljF,  POCIIC,  POF3 :  In  these  com¬ 
pounds  both  the  phosphorus  and  the  fluorine7  shielding 
increase  with  fluorine  substitution.  Pauling  has  pro¬ 
posed12  structures  for  the  phosphorus  oxyhalides  such  as 

0-  0 
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'  K.  Livingston,  J.  I’hv>.  Them.  57,  4%  (1953i. 

"  D.  \Y.  McCall  ami  H.  S.  Gntnuskv,  f.  Chem.  Phv>.  21,  1300 
( 1053). 


11  K.  S.  Pitzer,  J.  Am.  Chem.  Soc,  70,  2140  (1943). 

12  L.  Pauling,  The  Nature  0 1  the  Chemical  Bam .  (Cornell  Uni 
versify  Press,  Ithaca,  New  York,  194S),  second  edition,  p.  S3. 
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However,  in  order  to  give  the  observed  trend  in  the 
phosphorus  shifts,  there  must  be  double  bond  structures 
of  the  type 


o- 

P  and  perhaps 

/  \ 

Cl  F+ 

C1 

(2a) 


0 


Gi¬ 


ll 

P 

\ 

F+ 

C1 

(2b) 


The  increase  in  phosphorus  shielding  with  fluorine 
substitution  again  reflects  the  greater  ease  with  which 
fluorine  forms  double  bonds,  compared  to  chlorine. 
The  increase  in  the  fluorine  shielding  upon  continued 
substitution  is  analogous  to  that  observed  in  the  chloro- 
fluoromethanes  ;5  competition  among  the  fluorines 
makes  each  P— F  bond  less  double,  increasing  the 
fluorine  shielding;  but  in  the  polyfluoro  compound  the 
P—  F  bonds  still  have  more  double  bond  character  than 
P— Cl  bonds.  Although  the  effects  of  fluorine  substitu¬ 
tion  are  qualitatively  the  same  in  the  tri-  and  penta- 
valent  compounds,  the  magnitudes  are  considerably 
less  in  the  latter. 

(c)  P0C13,  PSCI3:  The  phosphorus  shielding  is 
higher  in  P0C13  than  in  PSCI3.  If  we  consider  the  an¬ 
alysis  above  as  valid,  this  provides  a  direct  verification 
that  oxygen  forms  double  bonds  more  readily  than  does 
sulfur  and  that  P  =  0  structures  are  important  in  the 
oxyhalides,  as  expected.11  By  analogy  with  the  trihalicles, 
if  the  ionic  forms  (la)  were  the  most  important  in  both 
POCl;iand  PSC!i,  the  greater  electronegativity  difference 
in  P—0  would  reduce  the  phosphorus  shielding  below 
that  in  the  sulfide.  The  double  bond  structure  (1  b)  must 
be  the  main  structure  in  POCR  to  account  for  the 
greater  shielding  observed  in  it. 

(d)  H2PO(OH),  HPO(OH)a,  PO(OH)»:  The  phos¬ 
phorus  shielding  increases  as  hydrogens  on  the  phos¬ 
phorus  are  replaced  by  OH  groups.  This  suggests  the 
importance  of  structures  like  P— 0  — H  and  P~=OH+. 
Confirmation  is  provided  by  the  shielding  of  the  protons 
attached  to  the  phosphorus,7  which  is  greater  in 
HPO(OII)-  than  in  H2PO(OH).  On  the  basis  of  the 
structures  postulated,  there  should  be  a  greater  nega¬ 
tive  charge  on  the  phosphorus  in  HPO(OH)»  than  in 
HjPO(OH),  reducing  the  ionic  character  of  the  P  — H 
bonds  and  increasing  the  proton  shielding  (i.e.  decreas¬ 
ing  the  proton  5  values'1),  as  observed. 

(e)  F3PO,  FsPO(OH),  FPO(OH)2,  PO(OH);1:  This 
set  of  compounds  is  similar  to  the  hydrogen-hydroxyl 


series  in  (d)  above  and  also  the  phosphorus  oxyfluoro- 
chlorides  in  (b).  The  phosphorus  resonance  could  be 
found  only  in  F2PO(OH)  and  PO(OH)3,  with  the  phos¬ 
phorus  shielding  greater  in  the  fluorine  derivative. 
Furthermore,  the  fluorine  shielding  increases  with 
fluorine  substitution.  These  trends  are  analogous  to 
those  in  ( b )  and  ( d )  and  are  explainable  in  the  same 
manner,  assuming  P_=  F+  electron  distributions  are 
relatively  more  important  than  P-=OH+. 

{f)  P4 ;  P4S3;  PH3:  As  remarked  in  the  initial  dis¬ 
cussion,  the  phosphorus  shifts  are  very  large  in  these 
compounds,  particularly  in  P4.  P4  and  P4S3  have  ab¬ 
normally  small  bond  angles;  the  associated  strain  no 
doubt  affects  the  bond  hybridization  and  could  be 
responsible  for  the  observed  effects.  The  high  value  in 
PH3  could  result  from  ionic  forms  P-H+,  which  increase 
the  charge  on  the  phosphorus  and  should  increase  its 
shielding.  A  sizeable  chemical  shift  exists  between  P4  (s) 
and  P4  (in  CS2)  indicating  there  are  significant  inter- 
molecular  effects  in  the  solid,  as  was  evidenced  earlier 
by  relaxation  of  the  infrared  selection  rules  in  the  solid.13 

In  conclusion  it  is  instructive  to  make  an  intercom¬ 
parison  of  the  influence  of  bond  type  upon  the  mag¬ 
netic  shielding  of  phosphorus,  fluorine,  and  hydrogen 
nuclei.  In  the  case  of  phosphorus  we  have  seen  that  an 
increase  in  double  bond  character  increases  the  nuclear 
shielding,  while  for  fluorine  the  shielding  is  decreased.3 
In  a  double  bond  the  charge  distribution  is  P~=F+; 
the  changes  in  net  formal  charge  are  opposite  for  the 
two  nuclei,  as  are  the  associated  chemical  shifts.  On 
the  other  hand  an  increase  in  ionic  character  decreases 
the  phosphorus  shielding  while  increasing  that  of 
fluorine ;  but  again  this  is  consistent  with  the  opposite 
charges  that  ionic  character  places  upon  the  two  nuclei, 
P+F-.  Hydrogen  does  not  form  double  bonds,  and  its 
shielding  decreases  with  an  increase  in  ionic  bond 
character,  R“H+;  however,  there  are  appreciable  bond 
hybridization  effects  as  well.8,8  The  chemical  shifts  in 
the  magnetic  shielding  of  all  three  nuclei  appear  there¬ 
fore  to  correspond  qualitatively  to  what  we  might  call 
differences  in  “electron  density”  about  the  nuclei. 
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SEVERAL  attempts  have  been  made  to  determine  the  spin  of 
Si59  with  either  inconclusive,1  conflicting,2  or  rather  indirect3 
results.  The  consensus  favors  the  value  1/2.  Nuclear  shell  theory 
gives  1/2,  3/2,  and  5/2  as  the  only  possibilities.1  We  have  elimi¬ 
nated  values  greater  than  1,  thereby  establishing  the  spin  of  Si2* 
as  1/2.  This  was  accomplished  by  resolving  the  multiplet  struc¬ 
ture  produced  in  the  fluorine  nuclear  magnetic  resonance  by  the 
electron  coupled  interaction4  of  the  Si22  spin  with  the  F19  spins, 
using  isotopically  enriched  SiMF,. 

In  pure  liquid  or  gaseous  Si29F,,  the  fluorine  resonance  should  be 
a  multiplet,  with  2/ (Si29)  + 1  equally  spaced  components  of  the 
same  intensity.4  The  SiF,  actually  used  was  prepared6  from  a 
sample  of  SiCh  whose  composition  was  given  as  30,  69,  and  1 
percent,  respectively,  of  Si28,  Si29,  and  Si30.  The  enriched  sample  of 
SiC>2  was  supplied  by  the  Oak  Ridge  National  Laboratory.  The 
spins  of  Si28  and  Si30  are  presumably  zero  so  the  Si28F,  and  Si^F, 
should  give  a  single  fluorine  resonance  with  31  percent  of  the  total 
absorption,  located  at  the  center  of  the  Si29F,  multiplet.  The  high- 
resolution  spectrometer  and  experimental  procedures  used  have 
been  described  elsewhere.0  A  typical  oscilloscope  photograph  of  the 
fluorine  resonance  observed  in  the  liquefied  sample  of  SiF,  is 
reproduced  in  Fig.  1. 

Figure  1  shows  a  symmetrical  multiple  resonance  of  three 
components  with  approximate  peak  intensities  of  1:2:1.  The 
separation  of  the  two  outer  components  is  46  milligauss.  The 
number  of  lines  observed  is  compatible  with  a  spin  of  either  1/2 


Fig.  1.  The  F19  nuclear  magnetic  resonance  in  a  mixture  of  Si!,F,  and 
Si"F,.  The  center  line  is  from  the  Sia8F,.  The  two  outer  lines  are  from  the 
SiMF,;  their  separation  is  46  milligauss.  A  sawtooth  sweep  of  135  milligauss 
was  used ;  the  pip  on  the  far  right  is  from  the  return  trace. 


or  1.  For  a  spin  of  1/2  the  predicted  integrated  absorption  of  the 
components  should  be  34.5:31:34.5,  and  for  a  spin  of  1, 23:54:23; 
in  the  latter  case  the  Si28F,  and  SiMF,  lines  overlap  the  center  com¬ 
ponent  of  the  Si29F<  triplet.  The  observed  peak  intensities  corres¬ 
pond  more  closely  to  the  spin  1  case  than  to  the  spin  1  /2.  However, 
a  spin  of  1  is  eliminated  by  the  shell  theory;  moreover,  there  is  a 
convenient  mechanism  to  explain  the  intensities  observed. 

In  aqueous,  concentrated  H2Si29F#  the  fluorine  resonance  was 
not  split  by  the  Si29  spin  but  was  broadened  to  about  4  milligauss 
compared  to  2  milligauss  in  H2Si28F,.  This  effect4  is  produced  most 
likely  by  fast  chemical  exchange7  involving  the  fluorines.  In  the 
case  of  SiF,,  hydrolysis  would  result  from  a  trace  of  water,  giving 
HF  which  would  catalyze  exchange  of  the  fluorines.  Depending 
upon  the  exchange  rate,  there  can  be  a  significant  reduction  in  the 
peak  intensities  and  some  broadening  of  the  component  lines.  In 
Fig.  1  noise  obscures  any  broadening  of  the  outer  components,  but 
this  explanation  is  supported  by  the  observation  that  the  ratio  of 
the  intensity  of  the  outer  components  to  that  of  the  center  line 
decreased  at  higher  temperatures,  for  which  chemical  exchange 
would  be  faster.  The  isotopic  composition  of  the  SiF,  was  checked 
by  an  infrared  analysis  which  gave  28±2  percent  Si28F,. 

Figure  1  suggests  an  amusing  and  direct  way  to  set  an  experi¬ 
mental  maximum  for  the  gyomagnetic  ratio  of  Si28.  The  splitting 
in  gauss  of  the  fluorine  resonance  in  Si'4F,  is  given  as4  SE  =  Cy 
X(Si A)/h,  where  C  is  a  constant  determined  by  the  electronic 
structure  of  SiF,.  The  central  component  in  Fig.  1  is  from  the 
Si28F,  and  splitting  of  it  greater  than  2  milligauss  would  be 
discernible.  If  we  take  the  ratio  of  the  equations  for  the  splitting 
in  Si28F,  and  Si"F,,  introduce  the  values  of  2  and  46  milligauss  for 
the  splitting  and  the  experimental  value1  for  7 (Si29),  we  find 
7(Si28)  <0.05.  This  value  could  be  improved  by  taking  greater 
pains  to  obtain  maximum  resolution  and  to  minimize  the  effects  of 
chemical  exchange.  Similar  experiments  are  feasible  in  several 
other  cases,  including  C12H,  — C13H,  and  H2S32— H2S33  mixtures. 

t  Assisted  by  the  U.  S.  Office  of  Naval  Research  and  a  Grant  in -Aid  from 
du  Pont  and  Company. 
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The  proton  magnetic  resonance  absorption  has  been  measured  in  NH<C1,  NILBr,  and  NH<I  crystal 
powders  from  — 195°C  to  room  temperature,  and  for  NH4C1  to  200°C.  Line  width  transitions  were  found 
at  —144,  —171,  and  about  — 198°C  in  the  chloride,  bromide,  and  iodide,  respectively.  The  N— H  distance 
in  the  ammonium  ion  was  determined  to  be  1.035±0.01A  from  the  second  moments  of  the  broad,  low- 
temperature  absorption  lines  observed  in  NHiCl  and  NH<Br,  with  allowance  for  broadening  by  modulation 
effects  and  narrowing  by  zero-point  torsional  oscillations.  The  relatively  narrow  room  temperature  line 
shapes  are  consistent  with  the  broadening  estimated  for  inter-NH<+  interactions,  the  intra-NHC  inter¬ 
actions  averaging  to  zero  over  the  hindered  rotational  motions  of  the  NH4+  ions.  No  appreciable  line  shape 
changes  were  found  at  the  X  temperatures.  An  electrostatic  calculation  is  made  of  the  potential  barriers  to 
rotation  of  the  NH<+  ion  giving  values  comparing  favorably  with  experimental  inferences.  The  results 
clarify  the  motions  which  narrow  the  absorption  line  at  low  temperatures  and  emphasize  the  dynamic  aspects 
of  the  order-disorder  process. 


1.  INTRODUCTION 

IN  papers  I1  and  II2  methods  were  presented  by  means 
of  which  nuclear  magnetic  resonance  absorption  line 
shapes,  widths,  and  second  moments  can  give  informa¬ 
tion  about  internuclear  distances  and  motions  in  solids. 
In  I  a  brief  analysis  was  made  of  the  broad  proton 
resonance  in  several  ammonium  halides  at  low  tempera¬ 
tures,  giving  a  value  of  1.025 A  for  the  N— H  distance  in 
the  ammonium  ion.  At  higher  temperatures,  the  proton 
resonance  lines  in  ammonium  salts  are  narrower,  sug¬ 
gesting3  low-frequency  rotational  motions  of  the  am¬ 
monium  ions.  This  article  includes  detailed  experi¬ 
mental  results  on  polycrystalline  ammonium  halides, 
at  temperatures  of  —  195°C  and  higher.  The  determina¬ 
tion  of  the  N— H  distance  in  the  ammonium  ion  is 
treated  more  exactly,  and  the  motions  responsible  for 
narrowing  the  resonance  lines  are  considered.  In  addi¬ 
tion,  the  article4  following  this  one  gives  results  ob¬ 
tained  with  a  single  crystal  of  ammonium  chloride,  for 
which  a  more  complete  analysis  and  more  comprehen¬ 
sive  conclusions  are  possible. 

Since  the  discovery  by  Simon5  in  1922  of  the  anomaly 
in  the  specific  heat  of  ammonium  chloride  at  —  30.4°C, 
an  enormous  amount  of  experimental  data  and  theo¬ 
retical  discussion  has  been  published  on  the  properties 
and  nature  of  this  and  similar  transitions  in  other  am¬ 
monium  salts.  The  infrared  spectra  of  several  of  the 
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ammonium  halides,6'7  as  well  as  the  analysis  by  Lawson8 
of  the  specific  heat,  indicate  strongly  that  the  room 
temperature  modifications  contain  NH4+  tetrahedra 
which  are  randomly  oriented  in  the  two  possible  equi¬ 
librium  positions  in  the  unit  cells,  with  the  hydrogen 
atoms  directed  essentially  toward  four  of  the  sur¬ 
rounding  eight  halide  ions  at  the  unit  cube  corners. 
In  the  low-temperature  phase,  the  NH4+  orientations 
are  ordered,  suggesting  that  the  transitions  are  of  an 
order-disorder  type.  Conclusive  confirmation  of  this 
model  has  since  been  provided  by  the  more  direct  neu¬ 
tron  diffraction  studies9'10  of  the  structures  of  the  hydro¬ 
gen  and  deutero-ammonium  chlorides  and  bromides. 
However,  the  details  of  the  order-disorder  process  are 
still  not  understood,  and  in  particular  interesting  ques¬ 
tions  remain  concerning  the  motion  of  the  ions. 

Prior  to  our  earlier  nuclear  resonance  studies  of  the 
ammonium  salts,  qualitative  observations  on  NH4C1 
were  made  by  Bitter11  and  by  Alpert.12  Since  then, 
Sachs  and  Turner13  have  measured  proton  spin-lattice 
relaxation  times  in  NH4C1,  NH4Br,  NH4I,  and  NH4IO3 
from  —195  to  25°C.  Soutif14  has  reported  line  width 
versus  temperature  data  on  a  single  crystal  of  NH4C1, 
and  Cooke  and  Drain16  have  made  some  line  width  and 
relaxation  time  observations  on  NH4CI,  NH4Br,  and 

6  E.  L.  Wagner  and  D.  F.  Hornig,  J.  Chem.  Phvs.  18,  296,  305 
(1950). 
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ln  G.  H.  Goldschmidt  and  G.  D.  Hurst,  Phys.  Rev.  83,  88 
(1951);  86,  797  (1952). 

11  Bitter,  Alpert,  Poss,  Lehr,  and  Lin,  Phvs.  Rev.  71,  738  (1947). 

11 N.  L.  Alpert,  Phys.  Rev.  75,  398  (1949). 

13  The  results  are  summarized  by  E.  M,  Purcell,  Physica  17,  282 
(1951);  details  are  given  by  A.  M.  Sachs,  Ph.D.  thesis,  Harvard 
University,  1949. 

14  M.  Soutif,  Rev.  sci.  89,  203  (1951). 

13  A.  H.  Cooke  and  L.  E.  Drain,  Proc.  Phvs.  Soc.  (London) 
A65,  894  (1952). 
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NH<I.  The  results  generally  agreed  upon  can  be  sum¬ 
marized  as  follows.  For  each  ammonium  salt  there  is  a 
characteristic  temperature  range  above  which  the 
resonance  absorption  line  is  relatively  narrow  and  struc¬ 
tureless  and  below  which  the  line  is  very  broad  with 
unresolved  fine  structure.  The  resonance  absorption  does 
not  change  significantly  at  the  X  temperatures.  The 
width  of  the  broad  proton  resonance  arises  from  the 
magnetic  dipole-dipole  interaction  of  a  given  proton 
with  the  other  three  protons  in  the  same  NH<+  ion  and 
also  with  all  other  neighboring  magnetic  nuclei  in  the 
crystal.  A  conclusion  supported  by  the  present  work  is 
that  the  narrowing  of  the  line  marks  the  disappearance 
of  the  intra-NH4+  interactions  and  a  slight  modification 
of  the  smaller,  inter-ionic  interactions.  The  disappear¬ 
ance  of  the  intra-NH<+  interactions  results  from  random 
reorientations  of  the  NH,t+  ions  at  an  average  frequency 
comparable  to  the  inverse,  proton  spin-spin  interaction 
time  (~10“  cps). 

This  frequency  is  much  too  low  to  contribute  ap¬ 
preciably  to  the  specific  heat  or  most  of  the  other 
physical  observables,  and  one  main  advantage  of  nu¬ 
clear  magnetic  resonance  is  that  the  presence  of  such 
low-frequency  processes  can  be  detected.  The  measure¬ 
ments  raise  the  problem  of  interpreting  the  low-fre¬ 
quency  motions  in  terms  of  the  crystalline  potential 


temperatures.  The  absorption  intensities  are  arbitrary.  The  pairs 
of  vertical  lines  represent  the  field  modulation. 
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GAUSS 

Fig.  2.  Derivatives  of  the  proton  magnetic  resonance  absorption 
lines  observed  in  the  ammonium  halide  powders  at  room  tempera¬ 
ture.  The  absorption  intensities  are  arbitrary.  The  pairs  of  vertical 
lines  represent  the  field  modulation. 

barriers.  The  barriers  are  calculated  electrostatically 
with  results  comparing  favorably  with  the  values  de¬ 
duced  from  the  observed  torsional  frequencies.6  The 
electrostatic  calculations,  relaxation  times,  and  line 
shape  studies  combine  to  give  a  physically  reasonable 
dynamical  model  for  the  rotational  motions  of  the 
NH<+  ions  and  for  the  order-disorder  process. 

2.  EXPERIMENTAL  METHOD 

The  data  reported  here  were  obtained  with  the  auto¬ 
matic  recording  rf  spectrometer  and  cryostat  described 
elsewhere.16  In  general,  these  observations  agree  with 
the  earlier  ones,1,3  the  main  differences  arising  from 
more  effective  temperature  control  in  the  present  ar¬ 
rangement.  A  fixed  radio-frequency  of  about  27.2  me 
was  used.  The  applied  magnetic  field  was  modulated  a 
small  fraction  of  the  line  width,  and  the  derivative  of 
the  proton  resonance  line  was  plotted  automatically  as 
a  function  of  applied  field  by  a  recording  potentiometer, 
preceeded  by  a  narrow-band  amplifier. 

The  samples  were  polycrystalline,  reagent  grade, 
from  commercial  sources.  An  arbor  press  was  used  to 
make  cylindrical  pellets  §  inch  in  diameter  and  £  inch 
long.  The  temperature  was  measured  with  a  copper- 
constantan  thermocouple  imbedded  in  a  small  hole 
drilled  in  one  end  of  the  pellet.  Complete  line  shapes 
were  obtained  at  fixed  temperatures  over  the  tempera¬ 
ture  range  investigated.  There  were  significant  differ¬ 
ences  in  the  line  shapes  observed  upon  warming  and 
cooling  through  the  line-width  transition  region.  The 
differences  could  arise  from  a  lag  in  attaining  tempera¬ 
ture  equilibrium,  most  likely  from  the  less  effective 
temperature  control  upon  cooling;  however,  there 
might  be  a  true  hysteresis.  In  any  event,  all  data  re- 
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Fig.  3.  Line  widths  of  the  proton  magnetic  resonance  absorp¬ 
tion  in  the  ammonium  halide  powders  plotted  against  tempera¬ 
ture.  The  curves  drawn  for  NHjCl  and  NHJir  were  calculated 
from  T i  results  as  described  in  Sec.  3c. 

ported  in  the  line-width  transition  region  were  obtained 
upon  warming. 

Second  moments  of  the  absorption  lines  were  obtained 
by  numerically  integrating  the  observed  first  deriva¬ 
tives.  The  errors  given  in  the  text  for  second  moments 
are  the  statistical  errors  of  the  observations.  In  addi¬ 
tion,  nonlinearity  and  systematic  error  in  the  biasing 
current-field  change  calibration  introduce  a  further 
uncertainty  of  no  more  than  one  percent  of  the  ob¬ 
served  second  moment  values.  Both  errors  are  taken  into 
consideration  in  calculating  the  N— H  bond  distance. 
The  experimental  errors  in  the  line  widths  depend  some¬ 
what  on  the  line  shape;  usually  the  statistical  and 
maximum  systematic  errors  together  were  no  more  than 
two  percent  of  the  values  reported. 

3.  RESULTS  AND  DISCUSSION 

The  absorption  line  shapes  observed  at  liquid  nitrogen 
temperature  are  given  in  Fig.  1  and  those  at  room  tem¬ 
perature  in  Fig.  2.  The  line  widths  5H,  defined  as  the 
separation  in  gauss  between  the  maximum  and  mini¬ 
mum  of  the  derivative  of  the  absorption  line,  are  plotted 
against  temperature  in  Fig.  3;  also  included  are  the 
various  phases  of  the  ammonium  salts  and  their  transi¬ 
tion  temperatures.  Second  moments  ASy,  obtained 
from  some  of  the  better  observed  derivatives,  are 
plotted  against  temperature  in  Fig.  4.  With  the  excep¬ 
tion  of  the  room  temperature  and  liquid  nitrogen  tem¬ 
perature  data,  the  points  in  Figs.  3  and  4  usually  repre¬ 
sent  single  observations.  The  line  widths  and  second 
moments  given  in  Table  I  for  those  two  temperatures 


are  averages  of  at  least  three  different  plotted  line 
shapes.  Ammonium  fluoride  is  not  included  in  this 
report  as  its  structure  differs  significantly  from  that  of 
the  other  halides. 

a.  Low  Temperature  Line  Shapes  and  the  N  — H 
Distance  in  the  Ammonium  Ion 

Figure  4  shows  that  the  second  moments  in  the 
chloride,  and  probably  the  bromide  as  well,  have  reached 
limiting  values  at  —  195°C.  These  values  no  doubt  corre¬ 
spond  to  “rigid  lattices”  and  can  be  used  to  evaluate 
the  N— H  distance  in  the  ammonium  ion.  Van  Vleck’s 
theoretical  result17  for  the  second  moment  may  be 
written  for  the  proton  resonance  in  the  ammonium 
halides  as1 

A/+2=  (9/32)g2/32{  E  (3  cos^y- 1)V„  "6 
/<»' 

+  (1/2)  E  (3  cos2#,-/,  1) V ,/+} 

i,k 

+  (1/12)27(7,+ 1  )g,v  E  (3  cos20y,-  DV+  (1) 

g  is  the  proton  gyromagnetic  ratio;  0  is  the  nuclear 
magneton;  the  indices  j  and  /  refer  to  the  four  protons 


Fig.  4.  Second  moments  of  the  proton  magnetic  resonance 
absorption  in  the  ammonium  halide  powders  plotted  against 
temperature. 

17  J.  H.  Van  Vleck,  Fhys.  Rev.  74,  1163  (<943). 
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Table  I.  Proton  absorption  line  widths  and  second  moments  observed  for  the  ammonium 
halide  powders  at  — 195  °C  and  25  °C. 


Compound 

Line  width,  SH 

25°C  -195°C 

Second  moment,  AHi* 

25  °C  -195°C 

N-H 

distance 

NH,C1 

4.8  gauss 

23.4  gauss 

3.35±0.05  gauss2 

49.5±0.5  gauss2 

1.038±0.004A 

NELfBr 

4.6 

23.7 

2.70±0.03 

48.0±0.5 

1.031±0.004 

NH4I 

3.8 

5.9 

1.55±0.02 

6.85±0.2 

in  a  given  NH4+,  k  to  all  other  protons  and  /  to  all  other 
magnetic  nuclei  in  the  sample ;  r3%  is  an  internuclear  dis¬ 
tance  and  9ju  is  the  angle  between  i3k  and  the  applied 
magnetic  field.  If  the  crystal  structure  is  known  except 
for  the  N— H  distance  r,  Eq.  (1)  is  essentially  a  one- 
parameter  equation,  and  r  can  be  computed  from  it  and 
the  observed  second  moment.  In  the  earlier  calculation1 
two  factors  were  neglected :  (a)  the  possible  broadening 
of  the  experimental  absorption  line  by  the  field  modula¬ 
tion;  and  (b)  the  effect  of  zero-point  torsional  vibrations 
of  the  ammonium  ion  in  the  calculation  of  the  theoretical 
second  moment.  It  turns  out  that  the  two  effects  nearly 
cancel,  giving  corrected  values  for  the  N— H  distance 
agreeing  with  the  original  results.1 

Relatively  large  modulation  amplitudes  were  used 
initially.  Perlman  and  Bloom18  showed  by  an  approxi¬ 
mate  method  that  this  introduces  an  appreciable  error 
in  the  second  moments.  Andrew19  has  given  a  direct 
calculation  of  the  effect,  obtaining  the  true  second 
moment  by  subtracting  /im2/4  from  the  observed,  while 
the  Perlman  and  Bloom  correction  is  hJ/3.  hm  is  the 
modulation  amplitude,  i.e.,  one-half  the  values  shown 
in  Figs.  1  and  2.  An  experimental  check  was  made  of 
the  correction  by  plotting  the  room  temperature  proton 
resonance  in  ammonium  bromide  as  a  function  of  modu¬ 
lation  amplitude,  at  amplitudes  less  than  the  line  width. 
The  resultant  second  moments,  which  are  averages  of 
several  plots  at  the  two  smallest  modulations,  are 
plotted  against  h,n  in  Fig.  5,  and,  for  comparison,  curves 
are  given  corresponding  to  hm2/ 4  and  hJ/3.  It  is  seen 
that  the  data  are  in  better  agreement  with  the  hm2/ 4 
correction  term.  It  is  of  interest  to  note  that  the  line 
width  83  was  4.6  gauss  independent  of  hm,  even  at  the 
maximum  2 hm  of  3.35  gauss.  The  second  moments  re¬ 
ported  here  and  in  the  following  article  as  well  have 
either  been  corrected  for  the  modulation  or  else  the 
modulation  amplitude  was  small  enough  that  the  cor¬ 
rection  is  well  within  the  other  experimental  errors. 

The  correction  for  the  zero-point  torsional  vibration 
is  obtained  by  averaging  the  angular  dependence  in 
Eq.  (1)  over  the  motion,  instead  of  assuming  a  com¬ 
pletely  rigid  lattice  and  averaging  only  over  the  random 
crystalline  orientations.  By  the  spherical  harmonic 


18  M.  M.  Perlman  and  M.  Bloom,  Phys.  Rev.  88,  1290  (1952). 
We  wish  to  thank  Mr.  Bloom  who,  prior  to  publication  of  their 
results,  pointed  out  to  us  the  importance  of  this  correction  in  our 
original  data. 

la  E.  R.  Andrew,  Phys.  Rev.  91, 425  (1953). 


addition  theorem,  we  can  write 
(3  COS20jk  l)torsional  motion 

=  (3  cos2<V- 1) <i (3  cos20—  1)),  (2) 

where  d3k°  is  the  equilibrium  value  of  0yt  and  0  is  the 
angle  between  r,-*  and  its  equilibrium  position.  0  is 
appreciable  only  for  the  intra-NH4+  interactions;  and 
as  the  ammonium  ion  behaves  like  an  isotropic  oscillator 
(see  Sec.  3d),  we  find 

3  h 

(3  cos20— 1)>=  (1-|02+  •  • .)  - 1- - -,  (3) 

8t rWV 


where  I  is  the  moment  of  inertia  of  the  ion  and  v  is  the 
torsional  frequency  in  cm-1.  The  correction  factor 
(1  —  e)  is  the  square  of  Eq.  (3)  and  introducing  it  into 
Eq.  (1),  averaging  over  the  crystal  powder  and  simpli¬ 
fying,  there  results 


A  tf22= 


5832 


8 

•|M — gv2 


81  920  15 


i SV-'U  — e) 


9  1 

+-/P2  £  fyr6+-//(17+ 1  )g&  £  ris-\  (4) 
80  i.k  15  i,f 


The  first  term  is  the  intra-NH.i+  broadening;  the  second 
gives  the  interactions  between  protons  in  neighboring 


Fig.  5.  The  second  moment  of  the  proton  magnetic  resonance 
absorption  in  NH,Br  powder  at  room  temperature  plotted  against 
modulation  amplitude  The  limiting  value  at  zero  modulation 
is  2.70  gauss1.  The  solid  and  dotted  lines  are  theoretical  values 
obtained  by  adding  A„J/4  and  //„V3,  respectively,  to  2.70. 
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NH4+,  and  the  last  represents  the  effects  of  the  halogen 
and  nitrogen  nuclei  external  to  the  given  NH4+. 

For  NH4CI,  the  second  terra  in  Eq.  (4)  was  found  to 
be  approximately  6.5  gauss2  by  calculating  correctly 
the  interactions  of  an  NH4+  with  its  thirty  two  nearest 
neighbors  and  replacing  the  remaining  ions  by  points 
on  a  simple  cubic  lattice.  The  NH4+  ions  were  assumed 
ordered  a-a-a-a;  an  ordering  a-b-a-b ;  where  b  is  the 
inverted  configuration  of  a,  would  give  values  about 
5  percent  higher.  The  results  depend  on  r,  which  was 
taken  to  be  1.04A,  and  on  the  size  of  the  CsCl-type  unit 
cell,  which  has  been  reported20  to  be  3.820A  at  —  195°C. 
The  last  term  in  Eq.  (4)  is  0.1  gauss2  for  NH4CI,  de¬ 
termined  by  computing  the  actual  distances  to  the 
nearest  eight  chlorine  and  six  nitrogen  nuclei,  replacing 
the  remaining  nuclei  by  points  on  CsCl  and  simple  cubic 
lattices,  respectively,  and  using  partial  summations 
readily  obtained  from  the  total  summations  given  else¬ 
where.21  The  torsional  oscillation  frequency  of  NH4+ 
has  been  reported6  to  be  390  cm-1  in  the  chloride.  With 
an  r  of  1.04A,  the  moment  of  inertia  of  NH4+  is  4.826 
X  10~'°  g  cm2.  These  values  in  Eq.  (3)  give  0.913  for 
(1  —  e).  As  a  final  result,  the  observed  second  moment 
of  49.5±0.5  gauss2  leads  to  a  value  of  1.038±0.004A  for 
the  N—H  distance. 

For  NE^Br,  the  calculation  of  the  inter-NH4+ 
broadening  is  complicated  by  the  tetragonal  structure9 
at  low  temperatures.  However,  the  distortion  from  a 
CsCl  structure  is  slight,  and  the  distances  are  such  that 
only  small  errors  are  introduced  by  assuming  the  CsCl 
structure,  for  which  the  arithmetic  is  much  simpler. 
A  CsCl  lattice  constant  of  4.1 1A  was  adopted,  the 
average  of  c+a/( 2)*  as  reported9  for  — 145  °C.  The 
ordering  of  the  NH4+  ions  appears9  to  be  such  that  they 
are  parallel  along  a  given  c  axis  of  the  tetragonal  struc¬ 
ture,  with  those  along  adjacent  axes  oppositely  oriented. 
This  corresponds  to  an  ordering  a-b-a-b  in  the  CsCl 
lattice  used.  The  large  lattice  constant  for  NH4Br 
permitted  the  proton-proton  distance  summations  to 
be  approximated  by  points  on  a  simple  cubic  lattice 
after  calculations  to  the  six  nearest  NH4+.  Otherwise 
the  computations  were  analogous  to  those  for  NH4CI. 
The  inter-NH4+  broadening  was  found22  to  be  3.9  gauss2 
and  the  broadening  by  bromine  and  nitrogen  nuclei 
outside  the  given  NH4+  0.4  gauss.2  The  NH4+  torsional 
frequency  is  reported6  as  319  cm-1,  giving  a  value  for 
(1  — <)  of  0.894.  Combining  these  results  with  the  ob¬ 
served  second  moment  of  48.0±0.5  gauss2  yields  a  value 
of  1.031±0.004  for  the  N—H  distance. 

The  1.038±0.004A  value  for  the  N—H  distance  in 
NH4CI  powder  is  within  experimental  error  of  the 
values  1.031±0.007A  and  1.033±0.004A  from  the 


*  L.  Vegard  and  S.  Hillesund,  Anhandl.  Norske  Videnskaps- 
Akad.  Oslo.  I.  Mat-Naturv.  Ki.  No.  8  (1942) ;  Chem.  Zentr.  1, 930 
(1943). 

51 H.  S.  Gutowskv  and  B.  R.  McGarvey,  J.  Phvs.  Chem.  20, 
1472  (1952). 

“These  values  are  somewhat  smaller  than  those  given  in  I, 
which  wore  based  on  a  smaller  lattice  constant  4.03A. 


(1,0,0)  and  (1,1,0)  single  crystal  orientations.4  The 
1.031±0.004A  value  from  NH4Br  powder  appears  to 
be  identical  with  those  from  NH4C1.  These  proton  mag¬ 
netic  resonance  values  are  to  be  compared  with  the 
neutron  diffraction  results9  of  1.03±0.02A  for  both 
ND4C1  and  NH4C1  and  for  ND4Br.  Bernal23  has  implied 
that  a  change  in  environment  of  the  NH4+  ions  might 
produce  a  change  in  bond  length.  However,  any  such 
effect  is  small,  even  in  NH4Br  where  the  tetragonal 
structure  suggests  the  importance  of  Br_-proton 
attractions. 

b.  Resonance  Absorption  at  Temperatures  Above 
the  Line-Width  Transitions 

The  relatively  narrow  absorption  lines  shown  in  Fig.  2 
persist  without  significant  change  from  the  line-width 
transitions  in  the  neighborhood  of  —  150°C  to  room  tem¬ 
perature  and  higher.  If  the  motions  of  the  NH4+  ions 
associated  with  the  line-width  transition  were  restricted 
to  reorientations  about  a  single  specific  axis,  such  as  a 
Ci  or  C3  symmetry  axis,  the  intra-NH4+  broadening 
would  still  be  an  appreciable  fraction2  of  the  “rigid- 
lattice”  value.  Moreover,  for  the  NH4C1  single  crystal4 
the  angular  dependence  and  magnitude  of  the  room 
temperature  second  moments  are  given  satisfactorily 
by  a  model  in  which  the  intra-NH4+  broadening  is 
averaged  to  zero  by  reorientations  effectively  about  all 
three  of  the  C2  axes.  The  observed  line  broadening  is 
mainly  from  inter-NH4+  interactions  which  are  reduced 
only  partially  by  the  reorientations,  with  a  small  con¬ 
tribution  from  proton-halogen  interactions.  This  model 
predicts  a  value  of  about  3.40  gauss2  for  NH4C1  powder, 
with  disordered  orientations  of  the  NH4+  ions ;  the  ob¬ 
served  second  moment  is  3.35±0.05  gauss2. 

A  similar  model  also  appears  reasonable  for  NH4Br 
and  NH4I.  The  distance  between  NH4+  ions  is  large 
compared  to  the  N—H  distance  so  the  numerical  values 
for  NH4CI  given  in  Eq.  (8)  and  Table  VIII  of  the  follow¬ 
ing  article4  should  yield  good  approximations  for  the 
other  halides  simply  by  using  their  unit  cell  dimensions. 
In  this  manner  the  inter-NH4+  broadening  from  NH4Br 
with  a  unit  cell9  of  4.060A  is  2.50  gauss2.  The  observed 
room  temperature  second  moment  is  2.69±0.05  gauss2 
which  compares  very  well  with  the  calculated  value 
when  consideration  is  given  to  the  proton-bromine 
broadening  of  0.4  gauss2  in  the  rigid  lattice.  The  case  of 
NH4I  is  somewhat  different  in  that  the  CsCl  phase  is 
stable  only  between  —17.6  and  — 41.6°C,  with  a  unit 
cell  of24  4.37A.  Using  the  same  procedure  as  for  NH4Br, 
an  inter-NH4+  broadening  of  1.58  gauss2  was  estimated. 
This  is  to  be  compared  with  a  second  moment  of  1.75 
±0.05  gauss2  observed  at  —  35  °C,  with  allowance  for 
proton-iodine  broadening  which  is  about  0.5  gauss2 
in  the  rigid  lattice. 

Once  the  intra-NH4+  broadening  disappears,  the 

ni.  J.  M.  Bernal,  Proc.  Phys.  Soc.  (London)  A66,  S14  (1953). 

54  A.  Smits  and  D.  Tollenaar,  Z.  physik.  Chem.  B52,  222  (1942). 
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theoretical  second  moments  for  the  crystal  powders  are 
insensitive  to  the  state  of  motion  or  ordering  of  the 
NH.(+  ions;  and  even  in  a  single  crystal  the  angular 
dependence  of  the  line  shape  is  only  slightly  different 
for  spherical  rotation,  and  for  disordered  and  ordered 
reorientations.4  In  agreement  with  this  prediction,  the 
second  moments  observed  for  NH4CI  powder  are  con¬ 
stant26  between  3.35  and  3.40  gauss2  from  room  tem¬ 
perature  to  above  150°C.  However,  NH4CI  undergoes 
a  phase  transition  from  a  CsCl  to  a  NaCl  lattice  at 
184.3°C,  and  about  this  temperature  the  second  mo¬ 
ment  is  definitely  smaller,  for  instance,  2.75  gauss2  at 
186. 5°C.  Detailed  calculations  were  not  made  for  the 
NaCl  lattice,  but  the  observed  decrease  is  qualitatively 
consistent  with  the  change  in  lattice  structure  and 
dimensions.  In  the  case  of  NH4I,  the  transition  between 
the  CsCl  and  NaCl  structures  occurs  at  —  17.6°C. 
The  second  moment  in  the  NaCl  phase  at  room  tem¬ 
perature  is  1.55  gauss2  which  is  also  smaller  than  the 
value  found  in  the  CsCl  phase,  1.75  gauss2  at  — 35°C. 
Cooke  and  Drain16  have  mentioned  a  line  width  change 
from  4.5  to  5.8  gauss  at  —  17.6°C,  which  we  did  not 
observe  although  the  second  moment  changed  as  indi¬ 
cated.  However,  there  is  some  possibility  that  the  NaCl 
phase  supercooled  in  our  experiments. 

From  these  results  we  conclude  that  the  intra-NHC 
interactions  become  averaged  out  during  the  line- 
width  transition  and  remain  so  at  higher  temperatures. 
This  requires  that  the  reorientations  of  the  NH4+  ions, 
at  temperatures  above  the  line-width  transitions,  do 
not  remain  correlated  with  any  specific  axis  of  rotation 
for  times  greater  than  about  10~4  sec.  Accordingly, 
some  modification  is  required  in  the  model  postulating26 
one  dimensional  rotation  of  the  NH4+  ions  in  the  NaCl 
phase  of  NH4I,  on  the  basis  of  the  infrared  spectrum. 
Rotation  of  the  NH4+  ion  can  occur  about  one  of  its 
C 3  axes,  but  there  must  be  interchange  in  the  orienta¬ 
tion  of  the  axis  about  which  the  rotation  occurs.  The 
relatively  low  interchange  frequency  required  would 
not  affect  the  infrared  spectrum  or  heat  capacity.  Three 
models  for  the  NaCl  phase  of  ND4Br  and  ND4I  were 


Table  II.  Frequency  factors  and  activation  energies 
for  reorientations  of  the  NH«+  ion. 


('nM)HHinil 

T\  measurements* 

( rAo  E 

Line-wMtli  transitions 
TV'  Ii 

\H,ri 

2.  OX  10  11  sec 

4740  cal 

-i44.rc 

4740'  cal 

XHjKr 

2.0X10” 

32X0 

- 170.9 

3340 

NH,I 

i.3x  in- " 

2900 

-198.1'* 

2550 

*  Tllf*-e  \ahlf'  arefrtiiu  Sadi;',  reference  IS.  *4  —  l/2x( 7 An. 

h  At  this  temiierature  the  line  width  is  half  way  betwmi  the  limiting 
values. 

*  Assumed  in  calculating  h  values  fur  XH«Br  and  NHJ. 

,T  This  value  was  extrapolated  hy  os-aiming  the  line-width  transition  f««t 
NH4I  to  be  of  the  same  general  shape  as  f* *r  NtWl  and  NH^Br. 


26  In  addition  to  the  main  resonance,  with  a  width  of  4.8  gauss, 
several  samples  exhibited  a  narrow  central  line  gradually  increas¬ 
ing  in  intensity  at  temperatures  above  about  50'C;  the  sharp 
line  was  found  to  originate  in  adsorbed  moisture  which  could 
catalyze  narrowing  by  self-diffusion. 

-6  R  Plumb  and  I*  F.  Hornig,  J.  Them.  I’hys.  21, 3( 16  ( 1952 1 . 


found  satisfactory  by  Levy  and  Peterson9  in  their 
neutron  diffraction  studies;  of  these  the  single  approach 
model  appears  to  be  the  only  one  consistent  with  all 
the  available  evidence. 

c.  The  Line-Width  Transitions 

The  narrowing  of  the  absorption  line  depends  upon 
the  correlation  frequency  vc,  the  average  reorientation 
rate  of  the  NH4+  ions.  In  the  case  of  a  change  from  a 
broad  line  to  a  narrower  line,2  the  frequency  width  of 
the  line  5v  is  given  as  a  function  of  vc  by  the  equation 

(5v)2=B2+C2( 2/ir)  tan-1[a(5>'/Vc)];  (5) 

where  B  is  the  width  of  the  narrow  line,  (B-JrCr)i  is  the 
width  of  the  broad  line,  and  a  is  a  parameter  depending 
upon  the  line  shape  with  the  value  unity  for  a  Gaussian 
shape.  The  reorientations  are  presumably  thermally 
activated  so  the  variation  with  temperature  of  vc  is 
of  the  form 

Vc=ne-EIRTt  (6) 

which  when  combined  with  Eq.  (5)  should  describe  the 
line-width  transition.  In  principle  Eqs.  (5)  and  (6)  can 
be  combined  with  the  observed  line  widths  to  compute 
E,  the  activation  energy  required  for  reorientation,  and 
go,  the  frequency  factor.  However,  a  is  also  unknown, 
complicating  the  analysis,  so  a  simpler  approach  is 
followed. 

It  is  physically  reasonable  that  the  NH4+  reorienta¬ 
tions  should  be  the  most  effective  spin-lattice  relaxation 
process.  Sachs13  and  also  Cooke  and  Drain16  have  ob¬ 
served  the  temperature  dependence  of  the  spin-lattice 
relaxation  time  T j  and  thereby  evaluated  E  in  the 
different  ammonium  halides,  and  also  vo.  The  results  are 
summarized  in  Table  II.  These  E  and  vu  values  can  be 
used  to  compute  vc  values  in  Eq.  (6),  and  the  results, 
substituted  in  Eq.  (5),  should  reproduce  the  line-width 
transition.  The  line  shape  parameter  a  was  evaluated  by 
fitting  Eq.  (6)  to  the  T,  Sv,  and  vc  values  at  the  center 
of  the  transition,  giving  4.88  and  10.4  for  NH^l  and 
NH4Br,  respectively.  In  this  manner,  the  lines  drawn  in 
Fig.  3  were  computed  for  these  compounds.  The  fit  is 
generally  quite  good;  the  deviations  for  NH<Br  at  the 
ends  of  the  transition  are  outside  experimental  error. 
However,  the  line  shapes  change  during  the  transition, 
(compare  Figs.  1  and  2),  and  this  may  change  a  from 
the  assumed  constant  value.  In  view  of  this  uncer¬ 
tainty  the  agreement  is  as  good  as  one  might  expect 
and  is  surely  adequate  to  demonstrate  that  the  NH<+ 
reorientations  determine  T i,  as  well  as  the  line  widths  in 
the  transition  region. 

Wert  and  Marx27  have  described  a  procedure  which 
can  be  adapted  to  give  activation  energies  for  reorienta¬ 
tion  from  the  line-width  transition  temperatures  Tt. 
The  correlation  frequencies  vc  should  be  very  nearly  the 
same  at  the  line  width  transition  in  the  different 


27  C.  Wert  and  J.  Marx,  Acta  Metallurgica  1,  113  (1953). 
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halides,  and  also  the  frequency  factors  vD  are  only 
slightly  different.  Accordingly,  Eq.  (6)  may  be  written 
as  E—ATt  where  the  constant  A  equals  R  In  (Vo,/ v)  and 
can  be  evaluated  from  a  known  E  and  T The  E  value 
for  NH4CI  of  4740  cal  mole-1,  from  T ,  data,  combined 
with  the  observed  Tt  of  129°K,  gives  a  value  for  A  of 
36.8  cal  deg-1.  On  this  basis  the  transition  temperatures 
of  102.5  and  75°K  for  NH4Br  and  NH4I  give  activation 
energies  of  3340  and  2550  calories.  Small  corrections 
have  been  made  for  the  differences  in  v0,  using  the  Ti 
results;  if  the  infrared  torsional  frequencies  are  used 
for  vo,  the  results  are  not  affected  significantly.  The 
value  for  NH4I  from  the  transition  temperature  is  a  bit 
lower  than  the  2900  calories  obtained  by  Sachs  from  T\] 
however,  Cooke  and  Drain  suggest  from  their  T\  data 
an  approximate,  lower  value  of  1700  to  2200  calories. 

The  two  interpretations  given  above  depend  more  or 
less  upon  the  T\  results.  A  completely  independent 
evaluation  of  E  could  be  based  upon  the  second  mo¬ 
ments  observed  as  a  function  of  temperature,  which 
would  eliminate  the  uncertainties  associated  with  the 
line  shape  parameter  a.  This  is  touched  upon  indirectly 
in  a  discussion  given  in  the  following  paper.4 


d.  The  Motion  of  the  Ammonium  Ion 


The  experimental  evidence  presented  and  reviewed 
above  demonstrates  the  existence  of  hindered  rotational 
motion  of  the  ammonium  ions.  In  this  section  a  phe¬ 
nomenological  treatment  is  given  which  clarifies  the 
nature  of  the  rotational  motion. 

A  phenomenological  potential  for  the  ammonium  ion 
can  be  obtained  following  a  method  of  Nagamiya.28 
The  ammonium  ion  is  assumed  electrically  to  be  a 
system  of  four  point  charges  of  +e/4  at  the  vertices  of 
a  rigid,  regular  tetrahedron.  The  potential  energy  of 
each  proton,  at  a  point  in  a  CsCl-type  lattice, 

with  respect  to  all  charges  outside  of  its  own  ion,  is 


co+ct 


(7ry,+,  3\ 


where  d  is  the  unit  cell  size  and  r  is  the  N—  H  distance, 
taken  to  be  1.04A.  The  first  term  is  very  much  larger 
than  the  second,  but  the  second  has  an  angular  de¬ 
pendence  which  is  the  chief  hindrance  to  rotation  of  the 
ammonium  ion.  Therefore,  a  suitable  potential  for  the 
rotation  is 


Vo  1  /Xi*+  y,*+Zi*  3 

r=— £  ( - 

2  i-i  V  r4  5 


For  the  specific  model  of  the  ion  proposed  above 


IV 


35/e\2 
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(■>> 


!*  T.  Nagamiya.  Prnc.  Phvs.  Math.  Snc.  Japan  24,  1.17  110521. 


Table  III.  Potential  barriers  to  rotation  of  the  NH4+  ion  in  the 
ammonium  halides  (CsCl  modification),  evaluated  by  an  electro¬ 
static  calculation. 


Compound 

Vo  Uiq.  9) 

Vo  (ICq.  14 )» 

nh4ci 

390  cm-1 

5400  cal 

5250  cal 

6410  cal 

NH<Br 

319 

3730 

4180 

4645 

NH*! 

279 

2910 

2850 

4100 

*  Tiiese  values  are  obtained  from  the  observed  torsional  frequencies 
while  the  Vo  (Eq.  9)  results  are  entirely  from  an  electrostatic  calculation, 
b  Eo  =E  (from  7'i  in  Table  II)  +(3/2)/n-„b.. 


where  r,-  is  the  distance  from  a  given  nitrogen  nucleus 
to  a  chloride  ion  or  another  nitrogen  nucleus,  the  sign 
being  +  or  —  respectively.  Values  of  Vo  computed 
with  Eq.  (9)  for  the  different  halides  are  given  in 
Table  III. 

The  potential  in  Eq.  (8)  is  more  general  than  the 
particular  model  used  in  obtaining  Eq.  (9).  To  demon¬ 
strate  this  generality  the  summation  in  Eq.  (8)  is 
carried  out  by  introducing  new  coordinates.  An  x',  y',  s' 
system  is  fixed  in  the  ion  such  that  at  equilibrium  the  s' 
axis  coincides  with  the  s  cubic  axis,  but  the  x\  y'  axes 
are  rotated  by  x/4  with  respect  to  the  x,  y  cubic  axes. 
0,  <j>  are  the  polar  coordinates  of  the  s'  axis  and  ^  meas¬ 
ures  rotation  of  the  ion  about  the  s'  axis.  In  terms  of 
these  angular  coordinates  Eq.  (8)  becomes 

F=|V0{  -  (49/8)  cos40+  (21/4)  cos-0-  (21/40) 

+  (7/8)  sin40  (cos4i/'— cos4<£) 

+  (cos20+|  sin40)  cos4 <j>  cos4f 

—  |  (cos0+cos30)  sin4</)  sin+A).  (10) 

This  is  the  simplest  nontrivial  potential  having  the 
cubic  symmetry  of  the  lattice  and  the  tetrahedral  sym¬ 
metry  of  the  ion.  The  general  solution  of  the  Schrodinger 
equation  for  a  spherical  top  moving  in  such  a  potential 
is  formidable  and  we  confine  ourselves  to  the  limit  of 
small  oscillations.29 

For  small,  commuting  rotations  we  can  choose  as 
coordinates 


£=0  cos<£ 

i;=0sin0  (11) 

S'=<H-vH-x/4, 

the  angles  of  rotation  about  the  cubic  x,  y,  z  axes.  De¬ 
veloping  the  potential  energy  in  a  power  series  of  small 
rotations  about  the  equilibrium  position,  we  obtain 

Uni  12  32 

F= — | - +8(^H2+s2) - (£HV+.<-4) 

2  15  3 


-8(£y+*v+rm--- 


(12) 


Using  first-order  perturbation  theoiy  the  energy'  levels 


n  In  so  doing  we  give  up  the  possibility  of  calculating  tunnelling 
frequencies. 
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of  this  anharmonic  isotropic  oscillator  are  computed 
to  be 

6 

jE»1H2«3= - Pod-  (Kl“)-W'V-f'W3-f-3/2)/fW 


1  8 

- [E  (2«i2+3«j+5/4) 

16  i”i 

+  «l»2+»lM3+«2«3]  {fluY/V  0,  (13) 


where  hu=  (16Fo(A2/27))J  and  7  is  the  moment  of 
inertia  of  the  ammonium  ion.  If  v  is  the  observed  tor¬ 
sional  frequency,  i.e.,  hv—EiM—Em,  then  from  the 
energy  levels  we  deduce  that 


1  (hv+Sh2/2iy~ 

F„= - ■ 

16  /i2/27 


(14) 


The  experimental  torsional  frequencies,  from  spectral 
data,  are  listed  in  Table  III  with  the  corresponding  Fo 
values  calculated  using  Eq.  (14).  These  experiment 
based  Vo  values  agree  well  with  those  determined  by 
Eq.  (9).  So  the  simple  model  used  in  the  electrostatic 
calculation  of  Eq.  (9)  gives  reasonable  results,  support¬ 
ing  Pauling’s  view30  that  the  charge  on  a  polyatomic  ion 
is  distributed  on  the  surface  of  the  ion. 

We  turn  now  to  the  essential  meaning  of  Fo.  If  an 
NH4+  ion  rotates  from  one  equilibrium  configuration  to 
another  around  a  fourfold  cubic  axis  it  passes  over  a 
barrier  of  total  height  Vo .  On  the  other  hand,  if  an  ion 
rotates  around  a  threefold  cubic  axis  (with  one  proton 
fixed)  the  barrier  is  (16/9)  Fn,  as  may  be  seen  from 
Eq.  (8).  The  activation  energy  E ,  as  given  by  the  T\ 
and  line-width  transition  analyses,  corresponds  to  the 
energy,  above  the  ground  state,  of  the  NH4+  ions  in 


their  reorientations.  If  the  values  of  E  are  corrected  by 
adding  the  zero  point  torsional  energy  of  the  ground 
state,  the  total  barrier  height  Eq  is  obtained.  These 
experimentally  inferred  values  are  listed  in  Table  III; 
they  are  larger  than  Fo  but  less  than  (16/9)  Fo  and  thus 
are  between  the  barriers  for  rotation  about  a  cubic  four¬ 
fold  axis  and  a  cubic  threefold  axis.  Therefore,  the 
motion  of  the  ion  which  is  most  important  in  relaxing 
the  proton  spins  and  in  narrowing  the  absorption  line  is 
a  motion  which  classically  is  partially  free  and  partially 
hindered.  This  type  of  motion  accounts  for  the  averaging 
completely  to  zero  of  the  intra-NH4+  contribution  to  the 
line  width.  Moreover,  the  experimental  values  of  v» 
given  in  Table  II  are  the  right  order  of  magnitude  for 
classical  rotation  frequencies  of  the  ions. 

e.  The  Order-Disorder  Process 

The  preceding  results  and  conclusions  emphasize  the 
dynamic  nature  implicit  in  the  order-disorder  process 
at  T\,  as  was  pointed  out  some  time  ago  by  Purcell.13 
Even  at  low  temperatures  in  the  ordered  phase  the 
ammonium  ions  experience  occasional  random  reorienta¬ 
tions,  but  return  to  the  stable,  ordered  configuration. 
As  the  temperature  is  increased,  the  rate  of  the  re¬ 
orientations  and  the  fraction  of  ions  reorienting  at  a 
given  instant  increase.  The  slightly  lower  energy  of  the 
ordered  phase  results28  mainly  from  the  electrostatic 
interactions  of  neighboring  NH4+  ions.  When  too  large 
a  fraction  of  the  NH4+  ions  are  in  motion,  the  coherency 
of  these  electrostatic  interactions  and  the  return  of  the 
ions  to  the  ordered  configuration  break  down,  producing 
a  cooperative  order-disorder  process.  The  over-all  mo¬ 
tions  of  the  ions  are  not  modified  appreciably  by  the 
degree  of  order  primarily  because  a  reorienting  ion  is  in 
a  highly  excited  energy  level,  while  the  order-disorder 
energy  is  small,  i.e.,  RT *«£. 


30  L.  Pauling,  J.  chim.  phys.  45,  142  (1048). 
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The  nuclear  resonance  absorption  spectrum  has  been  calculated  theoretically  for  a  tetrahedral  4-spin 
system,  with  I  of  i,  for  orientations  in  which  the  external  magnetic  field  coincides  with  a  two-  or  three-fold 
symmetry  axis  of  the  tetrahedron  or  with  a  tetrahedral  edge.  Comparison  of  the  theoretical  results  with  the 
proton  resonance  observed  at  — 195  °C  in  a  single  crystal  of  ammonium  chloride  demonstrates  that  the 
NH,+  ion  is  oriented  with  the  N— H  bonds  pointed  essentially  towards  corners  of  the  unit  cube.  An  N— H 
distance  of  1.032±0.005A  is  determined  from  the  low-temperature  second  moments  of  the  resonance  line 
with  the  applied  magnetic  field  in  the  (1,0,0)  and  (1,1,0)  directions.  At  room  temperature,  the  anisotropy  of 
the  resonance  results  from  inter  —  NH<+  interactions,  with  intra  —  NH4+  interactions  averaged  to  zero  by 
hindered  rotational  motions;  free  rotation  of  the  ammonium  ions  does  not  occur.  A  small  decrease  with  tem¬ 
perature  of  the  second  moment  of  the  resonance  line  at  the  X-point  agrees  with  the  order-disorder  model  for 
the  transition.  The  motion  of  the  NH<+  ion  is  discussed. 


1.  INTRODUCTION 

THE  preceeding  article1  gives  the  results  of  some 
proton  magnetic  resonance  absorption  experi¬ 
ments  on  ammonium  halide  crystal  powders.  However, 
the  data  obtained  from  a  crystal  powder  consist  of  a 
single  line  shape  and  its  second  moment.  Accordingly, 
more  than  one  structural  parameter  can  seldom  be 
evaluated  from  the  powder  data,  in  this  case  the  N— H 
distance  in  the  ammonium  ion.  But,  if  a  single  crystal 
is  available  (about  1  cm3),  the  angular  dependence  of 
the  resonance  absorption  provides  additional  data 
which  permit  a  much  more  detailed  analysis  of  the 
structure  of  the  sample  and  of  the  nuclear  magnetic 
dipole  interactions  themselves.  This  is  exemplified  by 
Pake’s  study2  of  a  gypsum  single  crystal,  and  also  by 
some  recent  research  on  K2HgCl4  •  H20  and  K2SnCl4  •  H20 
single  crystals.3 

In  the  case  of  the  ammonium  halides,  the  experiments 
by  Soutif4  and  this  work  are  apparently  the  only  mag¬ 
netic  resonance  experiments  in  which  single  crystals 
were  used.  The  present  work,  believed  to  be  more  com¬ 
prehensive  than  that  of  Soutif,  leads  to  the  following 
main  results:  (a)  the  orientation  and  size  of  the  NH<+ 
ion  are  determined  with  results  in  complete  agreement 
with  those  obtained  from  neutron  diffraction;6  (b)  at 
temperatures  for  which  the  line  breadth  has  reached  its 
narrow  value,  the  remaining  breadth  is  due  exclusively 
to  interionic  broadening.  The  first  conclusion  is  based 
upon  a  theoretical  calculation  of  the  absorption  spec¬ 
trum  for  the  tetrahedral  4-proton  system  in  the  am- 

*  Assisted  by  the  U.  S.  Office  of  Naval  Research, 
t  Presented  in  part  at  the  Washington  Meeting  of  the  American 
Physical  Society,  Mav,  1952;  R.  Bersohn,  Phvs.  Rev.  87,  226(A) 
(1952). 

1  Gutowskv,  Pake,  and  Bersohn,  J.  Phys.  Chem.  22, 643  (1954). 
>  G.  E.  Pake,  J.  Chem.  Phys.  16,  327  (1948). 

1  Itoh,  Kusaka,  Yamagata,  Kiriyama,  and  Ibamoto,  J.  Phys. 
Soc.  Japan  8,  293  (1953). 

4  M.  Soutif,  Rev.  sci.  89,  203  (1951). 

‘H.  A.  Lew  and  S.  W.  Peterson,  Phvs.  Rev.  86,  766  (1952); 

J.  Am.  Chem.  Soc.  75,  1536  (1953). 


monium  ion  and  a  comparison  with  the  low  temperature 
line  shapes  and  second  moments.  The  second  conclusion 
is  obtained  by  using  Van  Vleck’s  procedure6  to  calculate 
the  second  moments  for  several  models  and  comparing 
the  results  with  the  anisotropy  observed  in  the  room 
temperature  resonance. 

Relatively  little  direct  information  is  obtained  on  the 
order-disorder  transition  itself.7  Firstly,  an  ion  in  a  dis¬ 
ordered  position  differs  from  an  ion  in  an  ordered  posi¬ 
tion  only  by  an  inversion;  a  “disordered”  ion  therefore 
has  the  same  discrete  spectrum  as  an  “ordered”  ion. 
Secondly,  the  small  interionic  interactions  are  not  very 
sensitive  to  the  degree  of  order. 

2.  EXPERIMENTAL  METHOD  AND  RESULTS 

The  general  experimental  procedure  was  identical 
with  that  described  or  referred  to  in  the  preceding 
article.  The  single  crystal  on  which  the  reported  meas¬ 
urements  were  made  was  rectangular  in  shape,  14  by  14 
by  7  mm ;  it  was  supplied  by  Dr.  H.  A.  Levy  and  Dr. 
S.  W.  Peterson.  Some  preliminary  experiments  were 
performed,  with  similar  results,  on  a  smaller  crystal 
very  kindly  supplied  by  Mr.  Meyer  Bloom.  The  crystal 
■was  oriented  visually  in  the  applied  magnetic  field ;  line 
shapes  were  then  plotted  at  several  small  angular  dis¬ 
placements  to  ascertain  the  orientation  giving  the 
maximum  or  minimum  line  width,  and  this  orientation 
was  undisburbed  in  the  record  runs.  The  crystal  shape 
caused  some  difficulty  in  orienting  the  applied  magnetic 
field  in  the  (1,1,1)  direction;  it  was  primarily  for  this 
reason  that  line  shapes  were  not  obtained  as  a  function 
of  temperature  for  this  orientation,  particularly  as  the 
other  two  orientations  gave  data  sufficient  for  the 
present  purposes. 

The  observed  derivatives  of  the  proton  absorption 
fine  are  given  in  Figs.  1  and  2  for  the  applied  magnetic 

*  J.  H.  Van  Vleck,  Phys.  Rev.  74, 1168  (1948). 

7  An  erroneous  statement  to  the  contrary  was  made  by  one  of  us 
(R.B.)  in  the  first  reference  (t)  cited. 
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Fig.  1.  Derivative  of  the  proton  magnetic  resonance  absorption 
line  in  a  single  crystal  of  ammonium  chloride  at  —  195°C  with  the 
magnetic  field  applied  in  the  (1,0,0)  direction.  The  dashed  line  is 
the  result  of  a  theoretical  calculation  described  in  the  text.  The 
vertical  pair  of  lines  represents  the  field  modulation. 

field  in  the  (1,0,0)  and  (1,1,0)  orientations,  respectively, 
at  —  195°C,  and  in  Fig.  3  for  tire  (1,0,0),  (1,1,0),  and 
(1,1,1)  orientations  at  room  temperature.  The  dashed 
lines  in  Figs.  1  and  2  are  theoretical  line  shapes;  the 
details  of  the  comparison  are  discussed  later.  The  ex¬ 
perimental  data  in  Figs.  1  to  3  have  been  replotted;  a 
direct  reproduction  of  an  original  recording  for  the 
(1,1,0)  orientation  at  —  1 60°C  has  been  given  in  the 
article8  describing  the  instrumentation.  The  line  widths, 
defined  as  the  separation  in  gauss  of  the  maximum  and 
minimum  of  the  derivative,  are  plotted  in  Fig.  4  as  a 
function  of  temperature  for  the  (1,0,0)  and  (1,1,0) 
orientations.  A  similar  plot  of  the  second  moments, 
computed  numerically  from  some  of  the  better  line 


Fig.  2.  Derivative  of  the  proton  magnetic  resonance  absorption 
line  in  a  single  crystal  of  ammonium  chloride  at  —  195°C  with  the 
magnetic  field  applied  in  the  (1,1,0)  direction.  The  dashed  line 
is  the  result  of  a  theoretical  calculation  described  in  the  text. 
The  vertical  pair  of  lines  represents  the  field  modulation. 


*Guto\vskv,  Merer,  and  McClure,  Rev.  Sd.  Instr.  21,  644 
(1953). 


shapes,  is  given  in  Fig.  5.  Table  I  gives  values  of  the 
line  widths  and  second  moments  observed  at  —  195°C, 
room  temperature,  and  just  above  and  below  the 
X  point,  for  various  orientations  of  the  crystal  with 
respect  to  the  applied  magnetic  field.  These  values  are 
averages  of  at  least  three  separate  observations;  the 
experimental  errors  are  discussed  in  the  preceeding 
article.1 

3.  THEORETICAL  COMPUTATION  OF  THE  SPECTRUM 
OF  THE  FOUR  PROTON  SYSTEM 

(a)  Introduction 

The  most  striking  confirmation  of  the  assumed  orien¬ 
tation1  of  the  ammonium  ions  comes  from  examination 
of  the  actual  line  shape.  A  discrete  spectrum  for  a 
system  of  four  protons  is  calculated.  The  smaller  inter- 
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Fig.  3.  Derivatives  of  the  proton  magnetic  resonance  absorption 
line  in  a  single  crystal  of  ammonium  chloride  at  room  temperature 
with  the  magnetic  field  applied  in  the  (1,0,0),  (1,1,0),  and  (1,1,1) 
directions.  The  vertical  pairs  of  lines  represent  the  field  modulation. 

actions  with  other  neighboring  nuclei  are  then  intro¬ 
duced  approximately  by  smearing  the  discrete  lines  so 
that  the  second  moment  agrees  with  experiment.  This 
procedure  is  in  the  spirit  of  Take's-  and  Andrew’s9  work 
on  the  two-  and  three-proton  systems.  Although  the 
resulting  line  has  structure,  unlike  the  two-proton 
system  the  structure  is  unresolved  and  has  qualitative 
rather  than  quantitative  significance;  one  cannot  use 
a  separation  between  peaks  as  a  measure  of  an  inter- 
nuclear  distance. 

The  magnetic  Hamiltonian  for  the  four-proton 
system  is 

•30=  -Z gm,i+ Z  Ai, (I.-  1,-UiJu),  (1) 

t**I  i<; 

where 

A  a  =  rd‘>,T3[  (3/2)  cosfy,— }], 

*  E.  R.  Andrew  and  R.  Kersohn,  J,  Chem,  Rhys,  18, 159  (1950). 
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and  da  is  the  angle  between  the  field  H  and  r.y,  the 
vector  connecting  dipoles  i  and  j;  g  is  the  proton  gyro- 
magnetic  ratio  and  /3  the  nuclear  magneton. 

Because  we  are  interested  only  in  transitions  in  which 
angular  momentum  components  change  by  ±h,  only 
that  part  of  the  dipole-dipole  interaction  has  been  in¬ 
cluded  which  commutes  with  the  Zeeman  energy.  If  M 
is  the  quantum  number  associated  with  the  z  component 
of  the  total  angular  momentum,  we  expect  two  singlets 
(M  =  ±2),  tw'o  quartics- (M  =±1),  and  a  sextic  (if = 0) . 
Actually  the  sextic  will  factor  into  two  cubics  if  the 
states  are  labelled  by  I ,  the  total  angular  momentum 
quantum  number.  For  an  arbitrary'  orientation  of  the 
external  magnetic  field  no  further  factoring  is  possible. 


(b)  Conservation  Conditions  for  Factoring  the 
Secular  Equation 

For  the  purpose  of  seeing  whether  further  factoring  is 
possible,  a  representation  is  introduced  which  diagonal¬ 
izes  the  squares  of  the  operators 

L-I1+I2,  S=  I3— }— 14,  and  I=L+S. 

It  is  appropriate  to  introduce  the  quantum  number  /, 
because  P  commutes  with  all  terms  of  the  form  I, -It. 
The  matrix  elements  for  the  general  case  are  most  easily' 
calculated  in  an  AfiM 3Mt  representation  which 
makes  necessary'  the  transformation  (given  in  Table  II) 
from  the  LSIM  basis  to  the  basis. 

Table  III  shows  the  matrix  elements  of  the  Hamil- 


Fio.  5.  The  second  mo¬ 
ment  of  the  proton  mag¬ 
netic  resonance  absorption 
line  in  a  single  crystal  of 
ammonium  chloride  as  a 
function  of  temperature 
with  the  magnetic  field 
applied  in  the  (1,0,0)  and 
(1,1,0)  directions. 
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Table  I.  Line  widths  and  second  moments  observed  for  the 
proton  magnetic  resonance  in  a  single  crystal  of  NH,C1  at  different 
orientations  and  temperatures. 


Line  width,  $H  Second  moment,  AH22 

Temperature  in  Kauss  in  SaussI 

°C  (1,0,0)  (1,1,0)  (1,1,1)  (1,0,0)  (1,1,0)  (1.1.1) 


25 

8.2 

5.05  3.55 

6.70  ±0.10 

2.63  ±0.03 

-28 

8.3 

5.10 

6.62  ±0.15 

2.64  ±0.05 

-32  to  -36 

8.05 

5.10 

6.35  ±0.15 

2.57  ±0.05 

-195 

19.5 

23.3 

36.5  ±0.7 

54.6  ±0.6 

tonian  in  the  LSIM  representation.  Now  a  conservation 
law  for  any  set  of  angular  momenta  will  correspond  to 
the  factoring  of  the  secular  equations  because  they  are 
cubic  and  quartic  by  virtue  of  the  nonconservation  of 
the  angular  momentum  operator  D,  S2,  and  I2. 

The  Hamiltonian  can  be  written  down  in  the  form  of 
an  operator  which  will  commute  with  the  square  of  a 
given  angular  momentum  and  an  operator  which  will 
not  commute.  Sufficient  conditions  for  factoring  the 
secular  equation  will  result  from  setting  the  latter  oper¬ 
ator  equal  to  zero.  For  example  we  can  write  the  inter¬ 
action  Hamiltonian  in  the  form 

3C=  A  i2(Ii  •  I2— 3IxJzz)-\-Aw  (l3-l4—3/3Z/42) 

i 13+^4 23) C (Il+ 12)  •  I3—  3(7l2+72z)73  J 
2  C*4  i4~i~  -*4  24)  C  (hd"  I2)  •  I4— 3(7i2+/22)/42).  (2) 

+  2  (^4 13“ Zf23)[(ll—  I2)  'Is- 3(/lz—  /22)fW] 

+  2  (^4 14 — ^4  24)C(Ii — L)  -h— 3(7i2—  722)7<J 

When  Ai3=^23  and  Au=Au,  the  Hamiltonian  will 
commute  with  L2\  and  from  Table  III  one  can  see  that 
the  cubics  will  factor  and  the  quartic  will  factor  into  a 
cubic  and  a  singlet.  There  are  five  other  equivalent 
conditions  corresponding  to  the  six  different  ways  of 
choosing  a  pair  of  angular  momenta.  If  the  tetrahedron 
is  regular,  these  conditions  correspond  to  having  the 
magnetic  filed  in  one  of  the  six  planes  of  symmetry. 

If  one  asks  that  L 2  and  S2  be  simultaneously  con¬ 
served,  we  obtain  as  necessary  conditions 

A 13=  .4  14  =^23  =^24,  (3) 


with  two  other  equivalent  conditions.  For  a  regular 
tetrahedron,  these  conditions  correspond  to  having  the 
field  along  one  of  tire  three  twofold  axes  or  along  one 
of  the  six  edges.  In  this  case  we  find 

3C = A  ia(Ii  •  I2— 37  i272*) 34  (I3  •  I4— SIzzIiA) 

A-A  13(11+12)  •  (I3+I4)— 3(7i2+722)(732-|-742).  (4) 

If  S*—  (I1+72+I3)2  is  to  be  conserved,  then  we  have 

A 12= A 13=  A  23  (5) 

Au=Au=Asi, 

with  three  other  equivalent  conditions  corresponding 
to  orientation  of  the  field  along  one  of  the  four  threefold 
symmetry  axes  of  a  regular  tetrahedron. 

On  inspecting  Table  III  one  sees  that  the  conserva¬ 
tion  conditions  in  Eq.  (3)  result  in  the  factoring  of  the 
secular  equation  into  quadratics  and  singlets.  Equa¬ 
tion  (5)  is  equally  effective  although  it  is  not  obvious 
in  the  LSIM  representation. 

(c)  Spectra  for  Three  Orientations  of  the 
Magnetic  Field 

(i)  C3  Orientation 

For  the  C3  orientation  the  proper  representation  is 
LSIM  where  S=L+I3.  Tomita10  analyzed  the  nuclear 
resonance  spectrum  of  methane  by  an  independent 
study  of  the  four-spin  system.  His  paper  contains  the 
wave  functions  and  Hamiltonian  matrix  appropriate 
for  the  C3  orientation.  Tables  IV  and  V  give  the  energy 
levels  and  relative  intensity  of  the  ten-line  spectrum; 
in  these  and  subsequent  tables,  R  is  the  length  of  the 
tetrahedral  edge. 

(ii)  Ci  and  Edge  Orientations 

At  some  expense  of  complexity  in  notation  the  C2 
and  edge  orientations  of  the  external  field  can  be 
treated  together  because  they  are  both  cases  in  which 
L 2  and  S2  are  conserved  [condition  of  Eq.  (3)].  The 


Table  II.  Transformation  from  LSIM  basis  to  basis.  a<  and  ft  are  spin 

functions  for  an  individual  nucleus  t. 


2\pi 

2h 

\/2yq 

\/2y[j 


V2f7 

2y£, 

Hut 

2v/37n 


LSI 


1 

1 

2 

1 

'aiatcxtai 

1 

1 

2 

1 

1 

1 

1 

Piaxtjcu 

1 

1 

1 

-1 

-1 

1 

1 

atPtam 

1 

0 

1 

0 

0 

-1 

1 

aiaj/3 

0 

1 

1 

-1 

1 

0 

0 

aiawift 

1 

1 

2 

1  1  1 

1 

1  1 

ftdwrjcn 

1 

1 

1 

1  0  0 

0 

0-1 

aididja. 

1 

0 

1 

0-1  -1 

1 

1  0 

ftaaS»a< 

u 

1 

1 

0  1  -1 

-1 

1  0 

ftajayft 

0 

0 

0 

0-1  1 

-1 

1  0 

aiftajft 

1 

1 

0 

-2  1  1 

1 

1  -2 

aiajfSjft 

K.  Tomita,  Phys.  Rev.  89,  429  (1953). 


PROTON  RESONANCE  IN  NH.Cl 


655 


Hamiltonian  in  Eq.  (4)  is  equivalent  to 

3C=Ki4i2-^13)a2-3Z,2) 

+h(A3i~A13)  (S*-3S*)+Uu(.P-W.  (6) 

By  means  of  the  matrix11  {LSIM  |  LM lSM s),  one 
can  very  conveniently  calculate  the  matrix  elements  and 
the  energy  levels  with  results  which  are  given  in  Table 
III  and  VI,  respectively.  The  spectrum  consists  of  seven 
lines  if  the  magnetic  field  is  in  the  C2  orientation  and 
twenty  lines  if  it  is  parallel  to  a  tetrahedral  edge.  The 
intensities  of  the  lines  are  given  in  Table  VII.  Note 
that  Ajk  is  zero  for  a  regular  tetrahedron. 

4.  ABSORPTION  IN  NH4C1  BELOW  THE 
LINE-WIDTH  TRANSITION 

(a)  Comparison  of  Observed  Line  Shapes  with 
Calculated  Spectra 

From  the  plot  in  Fig.  5  of  second  moment  versus 
temperature,  it  is  clear  that  at  —  195°C  the  line  has 
reached  its  broadest  possible  value.  At  this  temperature 
motions  other  than  zero-point  vibrations  are  negligible 
and  we  will  assume  that  the  ammonium  ion  is  a  rigid 
torsional  oscillator. 

The  observed  and  calculated  absorption  derivatives  at 
—  195°C  for  the  (1,0,0)  and  (1,1,0)  orientations  of  the 
magnetic  field  are  shown  in  Figs.  1  and  2.  The  calcu¬ 
lated  derivatives  were  obtained  by  replacing  each  com¬ 
ponent  of  the  absorption  line  with  the  derivative  of  a 
Gaussian  function  whose  width  corresponded  to  the 
broadening  computed  in  Eq.  (7)  for  interactions  outside 
the  four-proton  system.  In  order  to  obtain  agreement 
between  the  experimental  and  computed  line  shapes,  it 
is  necessary  to  assume  that  the  three  C2  axes  of  the 
NH4+  tetrahedron  are  parallel  to  the  cubic  axes. 

That  the  agreement  is  not  quite  so  good  for  the 
(1,1,0)  curve  as  the  (1,0,0)  curve  is  perhaps  to  be  ex¬ 
plained  by  the  smearing  approximation  used.  The 
(1,1,0)  curve  consists  of  20  separate  transitions  as  com¬ 
pared  to  the  7  transitions  of  the  (1,0,0)  curve.  If  the 
different  transitions  were  smeared  differently  much 
better  agreement  could  be  obtained.  It  does  not  seem 
worth  while  to  pursue  this  point  because  NH4C1  is  a 
poor  example  to  use  for  this  comparison.  In  most 


Table  III.  Matrix  elements  in  LSIM  representation. 


M= 2 

(1122|77|  1122)=  —  lj> 

M=  1  (quartic) 

diagonal  elements 

(1121 1//|  1121)  =  ip  (1011  \H\  1011)= 

(UlHB\llll)=i(Au4-A3t)-}p  (0111 1 // 1 0111)  =  — ^ 
off-diagonal  elements 
(1121  |ff  1 1111)=  — A34) 

(1121 1//|  1011)  =  jj-2(A  13 1  A  a— An—  A  2<)=  3(1  111  |  H\  1011) 

(1121|//|01U)=i^(A23+A24-.413-A»)=-3(llll|ff|0111) 
(101 1 1 B 1 01 1 1) = \{A  i,-M  2,~ At,- A  u) 

M  —  0 

1— 0,2  (cubic) 

(11201/71 1120)=|/> 

(0000 1 Z7 1 0000) =(1100 1 // 1 1 100)  =  0 

(1120|/7iOOOO)=^(A14+A23-.424-J413) 

(1120 1 B 1 1 100)  =  ~  (3/1 12+3.4  U~P) 

(0000|77|  1100)=^(.421-.4„) 

7=1  (cubic) 

<01 10 1 77 1 01 10)  =  A  !4 ;  (1010 1  £7 1 1010)  =  A  ,2 
(1 1 10 1 77 1 1 1 10)  =  1^  —  (.4 12 +.4  34) 

(1110|/7|1010)=^(A14-t-.l24-.l,3-.423) 

(1010 1 77 1 01 10)  =  |  (.4 14+ .4  23-  A .4  24) 

(1 1 10 1 77 !  01 10) = ™  (,1 ,3+.  l  ,4-  A  tr -  A  24) 

P  =  ^i<hAjh 


crystals  where  the  XH4  or  XF4  groups  are  further  apart 
the  components  of  the  line  would  be  partially  resolved. 

In  the  case  of  solid  methane  the  question  has  been 
raised12  as  to  whether  or  not  the  molecules  exist  in 
ortho  (1=2),  mela  (1=1),  and  para  (1=0)  states  of 
total  proton  spin.  The  nuclear  magnetic  resonance 


Table  IV.  Energy  levels  when  external  field  is  along  a  C3  axis. 


State  in  LSIM  representation 

Energy 

<11221 

2g0H 

7=  (i+\/7/4)*(l!21 1  +  (1- vWdill  | 

S3ff+a+v'7/2)g»j3*7r> 

77=  (1— V7/4)KU21|  —  (i+V7/4)i(l]ll  I 

C077+  (» — x/7  /  2  )f?0tR-' 

ami,  <oim 

gPB— 

<li20|,  <0i00|,  <1100! 

0 

<11101 

-gVPR-’ 

(11101,  <01101 

1 

11  E.  LT.  Condon  and  G.  H.  Shortley,  The  Theory  of  Atomic  Spectra  (Cambridge  University  Press,  London,  reprinted  1051) 
M  T.  Nagamiya,  Progr.  Theoret.  Phys.  (Japan)  6,  702  (1051). 
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Table  V.  Intensities  of  the  lines  in  the  C3  orientation. 


AE  *»  ±  {Em  —Em  -1 HK 

Relative 

intensity 

±  (5/4— y/7/2)  ±0.073sW-s 

0.0744 

±2  s2^R~3 

0.2S00 

±  (-  HV7/2  WPR-*  =  ±  1  .OKgWR-3 

0.2534 

Mi+\/T/2)gV^3  =±1.57 3g^Hr3 

0.3716 

±  (5/4+V7/2  kVR~3=  ±2.573g2/327t-3 

0.0506 

data  for  methane,10  and  also  for  solid  hydrogen,13  are 
consistent  with  the  existence  of  nuclear  spin  isomers. 
Now  the  line  shape  agreement  in  NH4CI  depends  on 
the  existence  of  matrix  elements  between  states  of 
different  total  proton  spin,  which  proves  that  individual 
ions  do  not  exist  in  spin  states  of  definite  permutation 
symmetry.  In  view  of  the  exclusion  principle  the  normal 
vibrational  modes  must  be  such  that  the  motions  of 
neighboring  ions  are  highly  correlated  in  contrast  to 
the  molecules  of  solid  methane  or  solid  hydrogen. 

(b)  Second  Moments  and  the  N—H  Bond  Distance 

With  the  assumption  that  the  ions  are  ordered  so  that 
the  N—H  bonds  point  in  the  (1,1,1),  (1,1,1),  (1,1,1), 
and  (1,1,1)  directions,  Van  Vleck’s  formula6  for  the 
second  moment  becomes 

3,645  gW 

A  - - -  — {1—1  (X^+X-ri+X,4) }  (1  -  €) 

32,768  r6 

•W 

A - {1—  (Xi4+X^+X:i4)}  (1—  e) 

3  r° 

g2fi2 

+ — {8.7+28.0(Xi4+X.4+X*4)}+n.l  gauss.2  (7) 
d* 

Here  r  is  the  N—H  bond  distance,  g'  is  the  nitrogen 
gyromagnetic  ratio,  d  is  the  unit  cell  size,  and  X<  is  the 
direction  cosine  of  the  external  field  with  respect  to  the 
»th  cubic  axis.  The  factor  (1  —  «)  takes  into  account1  the 
torsional  oscillation  of  the  ammonium  ion.  The  first  two 
terms  in  Eq.  (7)  give  the  intra  —  NH<+  interactions. 
The  third  term  represents  the  interactions  between 
neighboring  NH<+  ions ;  the  coefficients  were  computed1 
assuming  r/d— 1.040/3.820.  The  0.1  -gauss2  term  is  the 
broadening  from  the  chlorine  and  nitrogen  nuclei 
outside  a  given  NH<+  ion. 

The  experimental  second  moments  at  low  tempera¬ 
ture,  given  in  Table  I,  when  substituted  in  Eq.  (7) 
result  in  N  —  H  distances  of  1.031±0.007  and  1.033 
±0.004A  for  the  (1,0,0)  and  (1,1,0)  orientations,  re¬ 
spectively.  These  values  are  in  excellent  agreement 
with  the  neutron  diffraction  values6  of  l.Q3±0.02A. 
This  agreement,  as  well  as  the  agreement  of  the  values 
for  the  two  different  crystal  orientations,  confirms  the 
orientation  assumed  for  the  NHt+  ion  in  the  unit  cell. 

>»  F.  Reif  and  E.  M.  Purcell,  Phvs  Rev  91,  631  (1953). 


5.  ABSORPTION  IN  NH,C1  ABOVE  THE  LINE- 
WIDTH  TRANSITION  TEMPERATURE 

(a)  Room  Temperature  Results 

Above  —  125°C  the  line  width  is  anisotropic  but  very 
slowly  varying  up  to  room  temperature.  The  theory 
given  in  Sec.  3  is  no  longer  applicable  as  the  measured 
second  moments  are  consistent  with  the  assumption 
that  the  line  broadening  is  due  solely  to  interionic 
interactions. 

Even  these  interactions  are  affected  by  the  motion 
of  the  ions,  and  one  can  try  using  the  dependence  to 
investigate  the  motion.  The  second  moment  in  general 
has  the  form 

AF22=g2/32d-6{c1(X14+X24+X34)-c2}+0.1  gauss2,  (8) 

where  the  second  term  is  an  estimate  of  the  contribu¬ 
tions  of  the  nitrogen  and  chlorine  nuclei.  The  angular 
dependence  of  this  term  has  been  neglected. 

Three  models  were  considered.  In  model  a  each  pro¬ 
ton  is  distributed  with  equal  probability  over  the  four 
vertices  of  the  tetrahedron;  in  model  b  each  proton  is 
distributed  with  equal  probability  over  the  four  vertices 
of  the  tetrahedron  as  well  as  over  the  four  vertices  of 
the  inverted  tetrahedron.  In  model  c  each  proton  is 
distributed  uniformly  over  the  sphere  circumscribing 
the  tetrahedron. 

Table  VIII  contains  the  values  of  a  and  c2  calculated 
for  each  model.  The  observed  values  at  room  tempera¬ 
ture  of  ci  and  c2  were  obtained  by  varying  the  data, 
given  in  Table  I  for  the  different  crystal  orientations, 
within  the  limits  set  by  experimental  errors  so  as  to 
obtain  a  completely  consistent  set  of  values  of  Ci  and  c2. 


Table  VI.  Energy  levels  when  external  field  is  along 
an  edge  or  a  Ci  axis. 


State  in  LSIM  representation 

Energy 

(1122 | 

i&m 

;=.vn(ii2ii-y12(iiii| 

X0H + {  (.1  n+.-i  j4 )+ /; 

77  =  ^12(1121  l+.5n(ll  11 1 

cSH+J  (3  is+.I  3<)~lt 

(10111 

tW-'-T 

<01111 

gftfr-'-y 

(0000| 

0 

s 

rH 

T 

II 

5^(.4m+.4h) 

7T'=4t{<1120|-|-(1100|) 

V  * 

~lhCll2+Au) 

(0110| 

■43< 

(10101 

■4 12 

(11101 

-UiS+.l,i) 

If  .113=.  1.,.  N.1*=1  If  1 11=  —2.1 

‘  11  162+2(82)» 

•Vi,=o 

A-1{10(.4„-.4,,H-44m.4„>» 
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A  comparison  of  the  values  in  Table  VIII  makes  it 
clear  that  the  second  moment  is  not  very  sensitive  to 
the  nature  of  the  motion  taking  place.  However,  sub¬ 
ject  to  this  limitation  the  room  temperature  second 
moments  show  that  not  only  is  the  crystal  disordered 
but  that  the  frequency  with  which  an  ordered  ion  be¬ 
comes  a  disordered  ion  (or  vice  versa )  is  greater  than  the 
frequency  associated  with  the  line  width,  which  is 
about  106  cps. 

(b)  Results  in  the  Vicinity  of  the  3w  Point 

The  great  changes  in  the  crystal  in  the  vicinity  of  the 
X  point  are  accompanied  by  almost  negligible  changes  in 
the  absorption.  The  disordered  model  b  shows  a  larger 
line  width  than  the  ordered  model  a  because  in  the 
disordered  crystal  the  average  proton-proton  distance 
is  smaller.  However,  when  the  crystal  becomes  dis¬ 
ordered,  it  also  expands  anomalously14  and  the  second 
moment  is  sensitive  to  the  size  of  the  crystal  through 
the  factor  d6.  The  two  effects  almost  cancel  so  that  the 
line  width  increases  only  slightly  through  the  phase 
transition. 

However,  Eq.  (8)  can  be  applied  to  the  observed 
second  moments  given  in  Table  I  for  temperatures  just 
below  and  above  tire  X  point.  With  3.84  and  3.85A 
assumed  for  the  unit  cell  sizes,14  c\  and  c2  values  of 

30.5  and  5.3  are  found  below  the  phase  transition  and 

32.6  and  5.9  above.  The  experimental  errors  are  such 
that  these  results  are  not  conclusive  in  themselves,  but 
upon  comparison  with  the  theoretical  values  in  Table 
VIII  they  provide  an  interesting  confirmation  of  the 
order-disorder  process. 


Table  VII.  Intensities  of  the  lines  in  the  edge  or  Ci  orientations. 


AE~±  {Em  —Em  - 1  ~kPH) 

Relative  intensity 

±Ci(-4i2+Aj«)-rA] 

±l-l(Avz+A:l)+hl 

hV 

±[(i-^)(AirfA,,)+x] 

A -5.1s 

±J^H - Y"^(-ll2+.43<)+I4J 

A5.fi 

±[(i+¥){'4,2+-43,w'] 

A5.fi 

A5.fi 

±[(5/4)C.4I2+.l,<)+A] 

}5,fi 

±  [(5/4)  (.4 12+ .4  34)— A] 

15.fi 

±i-4n 

1 

±^.434 

1 

If  .4,,=,l.n=-irW-3,  h=WPR~s 

(Ca  orientation) 

If  .4,2=  -2A>t=^R~>,  A  =  R 

"*  (edge  orientation) 

14  L.  Vegard  and  S.  Hillcsund,  Chem.  Zentr.  1,  930  (1943); 
P.  Dinichert,  Hclv.  Phys.  Acta  15,  462  (1943). 


Table  VIII.  The  anisotropy  of  the  second  moment,  described 
by  the  parameters  in  Eq.  (8),  computed  for  three  structural 
models  and  compared  with  the  room  temperature  measurements. 


C\  ft 

Model  a  (ordered) 

30.97 

4.86 

Model  b  (disordered) 

33.08 

5.65 

Model  c  (spherical  rotation) 

34.56 

6.14 

Experimental  values 

33.25 

5.88 

6.  THE  LINE- WIDTH  TRANSITION  AND  THE 
MOTIONS  OF  THE  AMMONIUM  ION 

The  observed  variation  with  temperature  of  the  line 
widths  and  second  moments  for  the  different  crystal 
orientations  suggests  that  the  motions  of  the  NH4+ 
ions  are  essentially  isotropic.  There  is  no  significant 
dependence  upon  crystal  orientation  of  the  character 
of  the  line-width  transition. 

The  change  in  line  width  with  temperature  was  cal¬ 
culated  for  the  (1,0,0)  and  (1,1,0)  crystal  orientations 
using  the  same  procedure,  frequency  factor  (vo),  and 
activation  energy  (E)  as  for  the  crystal  powder.1  The 
line  shape  factor  a  was  4.53  for  the  (1,0,0)  and  3.98 
for  the  (1,1,0)  orientation.  The  curves  drawn  in  Fig.  4 
give  the  results  of  the  calculation.  The  agreement  with 
experiment  is  reasonably  good;  certainly  there  is  no 
significant  difference  for  the  two  crystal  orientations. 

The  second  moments  measured  along  the  (1,0,0)  and 
(1,1,0)  crystal  axes  are  plotted  in  Fig.  5  against  the 
temperature.  An  analysis  of  these  data  discloses  the 
general  nature  of  the  motions  responsible  for  narrowing 
the  absorption.  If  we  let  St  be  the  contribution  to  the 
second  moment  at  absolute  temperature  T  of  the  proton- 
proton  interactions  on  a  given  NH4+,  then  it  is  found 
that  Sr/ So  is  independent  of  crystal  orientation,  within 
experimental  error.  In  Fig.  5,  the  curve  drawn  through 
the  transition  for  the  (1,0,0)  orientation  was  calculated 
from  the  St/So  values  for  the  (1,1,0)  orientation.  St  is 
given  by 

The  independence  of  S-r/So  upon  orientation  can  be 
derived  by  considering  the  motion  of  an  NH<+  ion 
during  the  time  required  for  rf  absorption.  In  particular, 
assume  a  model  in  which  any  given  edge  of  the  am¬ 
monium  ion  tetrahedron  is  at  its  initial  position  a 
fraction  (1  —  5a)  of  the  absorption  period  and  a  frac¬ 
tion  a  at  each  of  the  other  five  possible  positions,  i.e., 

A0-=(l -5a)  A ,-/’+  a  { £  *<  t  A °*  ;  ■ -  A 

=  (1  —  6 u)A\j  (9) 

where  An,  defined  as  for  Eq.  (1),  is  understood  to  be 
averaged  over  the  motion  of  the  ion  and  A%-  is  the 
value  of  Aij  when  no  motion  occurs.  Also,  we  have* 

Ar=  (9/4)£.<y  .1  ,/=  (9/4)  (1  -  6a)2£l<;,4  V. 
and 

S„=  (9/4 AW- 

Therefore,  we  obtain  the  result,  which  is  independent  of 
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orientation,  that 

ST/So=  (1-  6a)2.  (10) 

An  equivalent  physical  description  of  this  model  is  the 
following.  A  given  ion  spends  a  certain  fraction  of  its 
time  in  a  nonrotating  low-energy  state.  During  the  re¬ 
maining  time  it  is  randomly  excited  to  states  with 
higher  energy  in  which  reorientations  occur.  These 
reorientations  must  produce  in  the  high-energy  state 
a  uniform  spatial  distribution  of  the  ion  among  the 
possible  orientations  in  order  to  give  the  isotropic  result 
of  Eq.  (10).  The  temperature  dependence  of  the  second 
moment  arises  from  the  larger  fraction  of  time  spent  in 
the  high-energy  states  at  higher  temperatures. 

A  crude  classical  picture  is  that  the  ion  undergoes  a 
Brownian  motion  from  one  equilibrium  configuration 


to  all  other  equilibrium  configurations  with  the  same 
probability.  We  infer  that  the  state  or  states  in  which 
line  narrowing  occurs  must  lie  above  most  of  the  poten¬ 
tial  barriers  of  the  system.  The  same  conclusion  was 
reached  in  the  preceding  paper1  by  a  comparison  of  the 
potential  barriers  deduced  from  the  torsional  frequency 
and  from  the  variation  of  line  width  and  relaxation  time 
with  temperature.  Finally  we  note  that  the  isotropy  of 
St/So  implies  that  only  one  correlation  time  re  is  im¬ 
portant  for  this  frequency  range. 
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KIP,  Kittel,  and  co-workers  have  observed  electron  spin 
resonance  in  sodium1  and  in  lithium3  metal  at  microwave 
frequencies  and  at  300  Mc/sec.3  Measurements  on  these  metals 
at  lower  frequencies  have  been  made  by  Carver,  Holcomb,  and 
Slichterd  We  wish  to  report  the  results  of  some  independent  ob¬ 
servations,  at  frequencies  near  25  Mc/sec,  of  the  absorption  line 
shapes  for  lithium  from  77°  to  300 °K  and  for  sodium  from  77° 
to  385 °K. 

The  apparatus  was  identical  to  that  employed  for  broad  nuclear 
magnetic  resonance  lines5  except  for  two  concentre  pairs  of 
Helmholtz  coils  which  provided  the  main  magnetic  field  (0  to 
25  gauss)  and  the  field  modulation.  The  system  was  calibrated 
by  means  of  the  electron  spin  resonance  in  a,  a-diphenyl  /3-picryl 
hydrazyl.  The  lithium  and  sodium  samples  were  dispersions  pre¬ 
pared  by  high-speed  stirring  of  the  molten  metal  in  mineral  oil; 
particle  sizes  ranged  from  1  to  5<V,  averaging  15/t. 

The  observed  derivatives  of  the  absorption  lines  were  recorded 
automatically  as  a  function  of  applied  magnetic  field  at  a  fixed 
radio-frequency.  They  are  Lorentzian  in  shape  but  distorted 
somewhat  on  the  low-field  side.  This  distortion  can  be  attributed 
in  part  to  the  use  of  a  resonance  frequency  less  than  the  total 
half-width  of  the  absorption.  Thus,  there  is  still  absorption  at 
zero  applied  field;  moreover,  absorption  occurs  from  each  of  the 
opposing  circularly  polarized  components  of  the  oscillating  rf 
field.  Becker  and  Eisner5  have  reported  similar  detailed  results 
for  the  narrower  hydrazyl  resonance  at  frequencies  less  than  3 
Mc/sec.  Feher  and  Griswold3  have  pointed  out  the  possible  im¬ 
portance  of  skin  effects  in  distorting  the  line  shape. 

In  lithium,  the  resonance  at  room  temperature  has  a  width  of 
4.4  gauss  between  the  maximum  and  minimum  of  the  derivative. 
Widths  of  5  gauss  have  been  reported  by  Feher  and  Griswold3  and 
by  Slichter.4  The  Lorentzian  line  shape  suggests  that  spin-lattice 
relaxation  processes  determine  the  width,  which  corresponds  to  a 
T i  of  7 .5X10"*  sec.  This  relaxation  time  is  shorter  by  a  factor  of 
103  than  Overhauser’s  theoretical  prediction.7  Feher  and  Griswold3 
have  mentioned  that  the  line  width  increases  at  lower  tempera¬ 
tures;  in  our  sample  at  77°K  the  width  was  4.5  gauss,  only  slightly 
greater,  if  at  all,  than  at  room  temperature. 

On  the  other  hand,  the  width  of  the  resonance  in  sodium  ex¬ 
hibits  a  pronounced,  reasonably  linear  dependence  upon  tempera¬ 
ture,  as  shown  in  Fig.  1.  The  width  decreases  from  9.7  gauss  at 
385°K  to  3.6  gauss  at  77°K.  The  slope  of  0.018  gauss  deg-1  gives 
a  width  of  about  2  gauss  when  extrapolated  to  0°K.  The  results 
were  independent  of  the  resonance  frequency  between  20  and  30 
Mc/sec.  Feher  and  Griswold3  reported  a  width  of  10  gauss  at  room 
temperature  and  noted  a  temperature  dependence  as  well.  The 
different  temperature  effects  for  lithium  and  sodium  present  an 


Fig.  1.  The  temperature  dependence  of  the  electron  spin  resonance  ab¬ 
sorption  in  sodium  rri'-tal  dispersed  in  mineral  oil.  The  line  widtli  is  defined 
as  the  separation  between  points  of  maximum  slope. 


unresolved  problem.  Overhauser’s  analysis7  predicts  an  increase 
in  T i  and  therefore  a  decreasing  line  width  with  decreasing  tem¬ 
perature.  The  experimental  results  suggest  that  7)  includes  a 
temperature-dependent  and  a  temperature-independent  term,  the 
former  increasing  with  nuclear  charge.  Holcomb  and  Norberg* 
have  found  that  nuclear  relaxation  in  the  metals  is  sensitive  to 
impurities  and  imperfections,  which  could  contribute  to  the  elec¬ 
tron  relaxation. 

Electron  spin  resonance  was  observed  at  77°K  in  a  sample  of 
high-surface  sodium  prepared  by  mixing  molten  sodium  with 
surface-active  alumina  (Alcoa  F-20;  160  meter3  g_1).  The  pro¬ 
portions  of  sodium  and  alumina  were  such  that  if  the  sodium  film 
was  uniform  it  was  10A  thick.  The  line  shape  observed  is  identical 
with  that  for  the  dispersion  in  mineral  oil,  confirming  that  the 
resonance  is  from  the  sodium  and  not  the  supporting  medium. 
Unsuccessful  attempts  were  made  to  detect  the  electron  resonance 
in  K,  Rb,  Cs,  Cu,  Be,  Mg,  Al,  Ta,  Nb,  and  V  at  30  Mc/sec  at 
room  temperature  and  in  K,  Rb,  Cs,  and  Mg  at  77°K. 
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IT  is  conventional  to  use  the  formula  NH4OH  to  represent  the 
molecular  state  of  NH3  in  aqueous  solution.  This  formula  is 
most  likely  an  oversimplification.  In  fact  there  is  considerable 
argument  against  the  existence  of  discrete  NH4OH  molecules. 
Davis1  has  commented  on  the  problem  recently,  emphasizing  the 
weak  noncovalent  character  of  the  NH3— H2O  interaction,  and 
concluding  that  the  system  is  better  described  by  the  equation 

NH3+H2Oi=4NH,++OH-.  (1) 

Van  Velden  and  Ketalaar5  arrived  earlier  at  much  the  same 
conclusions. 

The  experimental  evidence  is  conflicting.  Rao3  found  that 
aqueous  ammonia  exhibited  the  same  Raman  spectrum  as  liquid 
NH3,  suggesting  no  new  important  molecular  species  are  present 
in  solution.  But  the  infrared  spectra  of  the  solid  ammonia  hydrates 
show4  absorption  bands  found  neither  in  ice  nor  in  solid  NH3. 
Early  thermodynamic  arguments5  were  based  on  the  partition  of 
NH3  between  H20  and  immiscible  organic  liquids  and  indicated 
that  half  or  more  of  the  NH3  remained  free  in  the  aqueous  solu¬ 
tions.  However,  Khakham6  reinvestigated  the  problem  by  this 
method  and  found  only  5  to  10  percent  free  NH3. 

An  independent  method  is  offered  by  nuclear  magnetic  reso¬ 
nance.  Chemical  shifts  in  the  resonance  position  are  a  sensitive 
indication  of  differences  in  the  electronic  environment  of  nuclei 
and  have  proven  useful  in  studying  the  structure  of  other  solu¬ 
tions.7  Figure  1  gives  the  proton  resonance  position  observed7  in 
aqueous  NH3  as  a  function  of  concentration,  at  room  temperature 
and  a  field  of  4200  oersteds. 


Fig.  1.  The  position  of  the  proton  resonance  in  aqueous  ammonia  as  it 
depends  on  concentration,  at  room  temperature.  3  is  the  shift  in  resonance 
position  referred  to  pure  HjO  and  pi  is  the  stoichiometric  fraction  of  protons 
as  NH.. 


The  resonance  at  all  concentrations  was  a  single  narrow  line 
(<1  millioersted).  If  discrete  molecular  species  existed  with 
lifetimes  longer  than  about7  10-3  sec,  separate  resonance  lines 
would  have  been  found  for  each  species.  The  conclusion  that  there 
is  fast  proton  exchange  agrees  with  the  inferences  of  Eigen8  and 
with  the  unsuccessful  attempts9  to  observe  directly  the  exchange 
rate  between  D20  and  NH3. 

When  chemical  exchange  produces  a  single  resonance  line,  its 
position  is  a  number  average  of  the  resonance  positions  for  the 
nuclei  in  the  absence  of  that  exchange.  If  Eq.  (1)  represented 
the  true  molecular  state  of  the  solutions,  the  resonance  position 
in  the  solutions  S ,  should  be  given  as7 

(2) 

and  defining  5H:o=0  as  was  done  experimentally, 

5„=/>i5nHi.  (3) 

pi  is  the  fraction  of  the  protons  in  NH3,  and  5nh.  and  5h.o 
represent  the  proton  resonance  positions  of  pure  NH3  and  H20. 
The  small  concentrations  of  NH.+  and  OH-  are  completely 
negligible  in  writing  Eq.  (2). 

The  observed  S,  values  in  Fig.  1  are  plotted  against  the  stoichio¬ 
metric  pi.  The  dotted  line  represents  Eq.  (3)  using  the  5nh. 
value  from  the  liquid.  The  large  deviation  of  this  line  from  the 
observed  data  is  strong  evidence  that  the  structural  model  of  free 
NH3  and  H20  is  inadequate.  Our  data  show  there  are  significant 
interactions  between  NH3  and  H2O,  but  the  state  is  a  very 
dynamic  one.  Eigen8  has  arrived  at  a  similar  model  from  supersonic 
absorption  measurements. 

The  importance  of  hydrogen  bonding  in  the  solutions  is  shown 
by  the  systematic  increase  in  the  5  values  of  NH3  (g),10  NH)  (1) 
and  the  aqueous  solutions.  However,  whether  one  wishes  to 
describe  the  state  of  association  as  NH(OH  molecules  or  NHr  HjO 
complexes  seems  at  this  point  a  matter  primarily  of  personal 
preference  or  convention.  We  are  attempting  to  learn  more  about 
the  different  molecular  species  actually  present,  by  cooling  the 
solutions  to  decrease  the  exchange  rate  so  separate  resonance  lines 
can  be  resolved. 
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EFFECTS  OF  RING  SIZE  ON  ELECTRON 
DISTRIBUTION  IN  SATURATED  HETEROCYCLIC 
COMPOUNDS* 

Sir: 

Measurements  of  the  basicities  of  cyclic  imines 
toward  trimethylboron,2  of  cyclic  ethers  toward 
chloroform  and  methanol-d,M  and  of  cyclic  sulfides 
toward  boron  trifluoride,6  have  shown  that  the 
basicity  changes  with  ring  size  in  the  order:  4  > 
5  >  6  >  3-membered  rings.  Two  differing  inter¬ 
pretations  of  these  observations  have  been  sug¬ 
gested. 

An  interpretation  based  on  steric  factors  was 
proposed  by  Brown  and  Gerstein2  to  account  for  the 
dissociation  of  the  addition  compounds  of  the  cyclic 
imines  with  trimethylboron.  According  to  this 
view  the  observed  order  of  basicity  results  from  a 
combination  of  F-strain  which  is  more  pronounced 
for  the  larger  rings,  and  I-strain  which  is  most 
important  in  the  3-membered  ring.  These  strains 
are  considered  to  be  steric  interactions  which  occur 
because  of  the  association  between  donor  and 
acceptor  molecules. 

On  the  other  hand,  steric  factors  alone  do  not 
account  adequately  for  the  results  of  the  hydrogen 
bonding  studies  with  cyclic  ethers3'4  nor  are  they  a 
likely  explanation  for  the  interaction  of  cyclic 
sulfides  with  boron  trifluoride.6  It  was  suggested 
that  the  basicity  differences  were  due  rather  to 
differences  in  electron  availability  caused  inherently 
by  the  different  sizes  of  the  rings — that  is,  the  elec¬ 
tron  distribution  is  altered  by  ring  size. 

Direct  physical  evidence  that  the  electron  dis¬ 
tribution  does  depend  significantly  on  ring  size  has 
now  been  obtained  by  observing  the  chemical 
shifts  in  the  proton  magnetic  resonance6  of  the 
cyclic  imines,  ethers  and  sulfides.  These  shifts 
are  a  sensitive  measure  of  differences  in  the  elec¬ 
tronic  environment  of  nuclei.  The  5-values6  ob¬ 
served  for  the  cyclic  compounds  are  given  in  Table 
I,  and  it  is  seen  that  there  are  relatively  large 
variations  with  ring  size  for  the  hydrogens  in  both 
the  a  and  fl  CHj  groups. 

In  each  series  of  cyclic  compounds,  the  5-values 
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Table  I 

Proton  Magnetic  Resonance  Shifts  in  Saturated 
Heterocyclic  Compounds 

Compound  . - 4-Value" - . 


«-CH, 

P-CHt 

(ch2)2o 

— 0,26 

(ch2),o 

—  .06 

— 0.25 

(ch2)4o 

-  .17 

—  .34 

(ch2)6o 

-  .16 

-  .39 

(CH2)S 

-  .23 

(CH2),S 

-  .19 

-  .19 

(CH2)4S 

-  .29 

-  .29 

(CH2)6S 

-  .29 

-  .34 

(CH2)sNH 

-  .37 

(CH2),NH 

-  .23 

-  .16 

(CH2)4NH 

-  .19 

-  .30 

(ch2)5nh 

-  .25 

-  .36 

•  Tlie  assignments  of  lines  to  the  a-  and  /3-CHj  groups 
usually  could  be  made  from  their  relative  intensities.  There 
are  some  uncertainties  in  values  for  the  cyclic  imines  because 
the  N-H  proton  resonance  was  not  always  resolved  from  the 
CHj’s. 

of  the  3-membered  rings  are  consistently  among 
the  most  negative  while  the  a-CHj’s  in  the  4-mem- 
bered  rings  are  among  the  most  positive.  It  has 
been  established6  that  more  positive  5-values 
correspond  to  lower  electron  densities  about  the 
proton.  Also,  a  low  electron  density  on  the  CH2 
groups  must  be  compensated  for  by  a  high  electron 
density  on  the  hetero  atom.  So  we  conclude  that 
the  hetero  atom  has  a  low  electron  density  in  the  3- 
membered  ring  and  high  in  the  4,  which  agrees 
with  the  observed  order  of  basicity. 

It  is  also  noteworthy  that  the  changes  in  electron 
distribution  with  ring  size  appear  to  be  more  sys¬ 
tematic  at  the  0-CH2  group  than  at  the  a.  As  the 
ring  size  increases,  the  5-value  for  the  /3-CH2 
group  approaches  monotonically  the  value  —0.39 
observed  in  cyclohexane.6 

Further  study  of  variations  in  electron  distribu¬ 
tion  with  ring  size  is  planned  with  these  and  other 
systems.  More  detailed  data  and  interpretations 
will  be  presented. 
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